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JET DIFFUSION IN LIQUID OF GREATER DENSITY 


By G. Abraham! 


SYNOPSIS 


A study on the influence of the density difference between a jetand the sur- 
rounding fluid on the diffusion by a circular vertical submerged turbulent jet 
is presented. The results are summarized in Figs. 1(a) and 1(b), which show 
how the density and the velocity along the axis of the jet depend on the initial 
density difference between the jet and the surrounding fluid and on the velocity 
of the jet when it issues from a nozzle. The theoretical basis for the families 
of curves of Fig. 1 is described together with laboratory experiments, which 
confirmed the theoretical considerations. The experiments were performed 
with water jets. The density of the jet was 1,000 kg per cu m (fresh water), 
the density of the surrounding fluid (salt solution) varied from 1,020 to 1,050 
kg per cu m. Because of the range of densities covered by the experiments the 
results apply to the case of sewage disposal in the marine environment. 


INTRODUCTION 


For the case of a turbulent jet surrounded by homogeneous still water, the 
jegree of dilution depends on the quantity and velocity of flow, the diameter of 
he outlet-pipe, and the density difference between the fluid of the jet and the 
surrounding water. The relationship between these factors is described here 
‘or the case of a turbulent jet issuing vertically upward in surrounding water 


vith a specific weight greater than that of the jet. Theoretical information is 


ee ee ee ——————— 
_ Note.—Discussion open until November 1, 1960. To extend the closing date one month, 
. written request must be filed with the Executive Secretary, ASCE. This paper is part 
£ the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 86, No. HY 6, June, 1960. 

1 Hydr. Engrg. Lab., Wave Research Projects, Univ. of California, Berkeley, Calif. 
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presented together with results of laboratory tests. The range of density dif- 
ferences covered by the experiments is such that the results obtained apply tc 
the case of sewage disposal in the marine environment. 


THEORY 


In the case of homogeneous surrounding water withno currents, the velocity 
u and the concentration c atany point of the jet dependin general on the follow- 
ing variables: 


x = the distance from the nozzle of the outlet-pipe measured along the 
axis of the jet; 


y = the distance from the axis of the jet; 

d = the diameter of the outlet-pipe; 

Ug = the velocity of the fluid of the jet at x = 0; 
Po = the density of the fluid of the jet at x = 0; 


Ap = Pg - Po: the density difference between the fluid of the jetat x = 0 anc 
the surrounding fluid; and 


v = the kinematic viscosity of the fluid of the jet at x = 0. 


The concentration c is defined as: 


in which p is the density at any point of the jet and pg the density of the sur- 
rounding fluid. 

Dimensional analysis shows? that the velocity u and the concentration c de- 
pend on the other variables as: 


Mee g (eX 40 
Uo (5. d’ Pog oF R | C.e F 0 e we. le te) wee, eh pr ere (2 
and 

Or fi eee 

Co (5. d’ Psn Fr. R| = © ie (Ss) 2 Wa, ce eee eee (3 
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and R is a Reynolds number, 
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The Reynolds number gives the ratio of the inertia forces to the viscous 
forces. Since this paper is restricted to turbulent jets, the Reynolds numbe! 


is considered to be sufficiently large to have no influence on the functions o 
Eqs. 2 and 3. 


ae : : ; 
2 “Diffusers for Disposal of Sewage in Sea Water,” by A M Rawn, F. R. Bowerma 
and N. H. Brooks, Proceedings, ASCE, Vol. 86, No. SA 2, March, 1960. i 
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The Froude number gives the ratio of the inertiaforces to the gravitational 
lorces. If at any point of the jet the local Froude number basedon the local ve- 
locity and the local density-difference with the surrounding water is large, the 
influence of the gravitational forces may be neglected at that point. It is evi- 
Jent that the inertia forces are less important at greater distances from the 
1ozzle. Hence it is possible that near the nozzle, the Froude number is suf- 
iciently large to neglect the influence of the buoyancy while it may be that at 
greater distances from the nozzle the buoyancy has to be considered. 

According to F. H. Schmidt? the problem may be divided into the following 
jarts: 


1. The non-buoyancy case, defined by F =~. For this case the buoyancy is 
1egligible at any point of the jet. 

2. The buoyancy case, defined by F>0O. For this case the initial momen- 
um of the jet is negligible. 

3. The intermediate case, defined by a finite value of F. For this case the 
nitial momentum is predominant near the nozzle, the buoyancy at larger dis- 
ances from the nozzle. 


A discussion of the various parts of the problem follows. 

The Non-Buoyancy Case, F =~.—Studies regarding the non-buoyancy case 
vere conducted by M. L. Albertson, Y. B. Dai, R. A. Jensen and H. Rouse, 4 
N. Forstall and E. W. Gaylord,® J. O. Hinze and B. G. van der Hegge Zynen,® 
;. Corrsin and M. S. Uberoi,” and W. R. Keagy; A. E. Weller, F. A. Reed and 
N. T. Reid.8 

These studies yield the following results: 


“or the zone of flow establishment iS <6 or 7) 


or y> $- C3 x, 


: d : 
ory<5 - C3 x; 


3 “On the Diffusion of Heated Jets,” by F.H. Schmidt, Sartoyck ur Tellus, Vol. 9, No. 

, August, 1957, pp. 378-383. / 
“Diffusion of Submerged Jets,” by M. L. Albertson, Y. B. Dai, R. A. Jensen and H. 
louse, Transactions, ASCE, Vol. 115, 1950, pp. 639-664, 

5 “Momentum and Mass Transfer in a Submerged Water Jet,” by W. Forstall and E. 
J, Gaylord, Journal of Applied Mechanics, Vol. 22, No. 2, June, 1955, pp. 161-164. 

6 “Transfer of Heat and Matter in the Turbulent Mixing Zone of an Axially Symme- 
ical Jet,” by J. O. Hinze and B. G. van der Hegge Zynen, Seventh Internatl. Congress 
»r Applied Mechanics, London, 1948, pp. 286-299. 
_74}urther Experiments on the Flow and Heat Transfer in a Heated Turbulent Air 
et,” by S. Corrsin and M. S, Uberoi, N.A.C.A. Tech, Note N 1865, April, 1949. 

8 “Mixing in Inhomogeneous Gas Jets,” by W. R. Keagy, A. E. Weller, F. A. Reed and 
/,T. Reid, Batelle Memorial Inst., The Rand Corp., Santa Monica, Calif., February, 
949, 5 re 


iy 
. 


peacennns 


é 


4 June, 1960 HY@ 


wash? apis che! ee Sa ee (Tb 
fo 
for y< s - C4 x. 
For the zone of established flow & 76 or 7) 
u __-k(y/x/? 
Tse 6.0 ys > sats hene. e Ber eee (8 
and 
u d 
ro =O y eee ee (9 
for 5 7Cy; 
ne Hy Kk (y/x)? ac 
Gute e At OS coe «ee ae 
and 
c 
m d 
rats = Co = Oleh Se ei <emne stew tS (pie eelaee er ene ee (i 
fe) 
for ; > CQ; 


in which u,, is the velocity at the axis of the jet (y = 0), cy, is the cor 


centration at the axis of the jet, and , k, Cy, Cg, C3 and C4 are dimensionles 
constants. To obtain a smooth transition between the zones it isnecessary th 


OD ee: 
Cz = 3 C; ee eae, Se 6 ol ee (1: 
and 
ae Tee | 
Cy = DY Cy oe cee Re See enw soe enn were (1: 


The experimental values for the constants yu, k, C, and Cg, and the corresponc 
ing values of C3 and C4 are shown in Table 1. 


The experiments of Corrsin and Uberoi’ show that in pee range 1< s aa 
(the range covered by their experiments), k = 96 [1 +0. 19(° —- S-1\} -2 and w=0 
= constant. As a consequence the constants k and no nat vary much fc 
L< e < 1.025, that is, within the range of - values that are of interestin sey 


age eae problems in the ocean. 
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Ps 
Hencefor 1< = < 1.025 we may describe the experimental results by aver - 


fe) 
ie values for k and wand according to theory by average values for Cy and 
29: Since water jets are involved in the disposal of sewage, the following have 
een Selected: k = 77, = 0.8, Cy = 6.4, Co = 5.2, Cg = 0.078, and C4 = 0.096. 
[hese are the values obtained in the experiments of Forstall and Gaylord.® 


TABLE 1.—PREVIOUSLY PUBLISHED DATA FOR NON-BUOYANCY CASE 


.Ibertson et al. 

‘orstall and water jets 

aylord salt solutions 

eagy et al. mixture of 
gasses 

linze and mixture of 

an der Hegge gasses 

“ynen 


Yorrsin and 
Jberoi 


The Buoyancy Case (F—+0).—Studies on the buoyancy case were made by B. 
2. Morton, G. Taylor and J. S. Turner? and by Rouse, C.S. Yih and H. W. 


fumphreys.! 
The theoretical work of Morton, et al.,9 yielded the following results: 
1/3 
5 9 oat 
Um =§% a( 3 a a) x LS 5 it hnaetea ee (14) 
nd 
§ (Pg -P 2.1/3 2 
oo) oe ae TEE aaa (15) 
Pg Geo AELD 


Ps - Po 
n which @ is a dimensionless constant and Q = 7 d? uy einer 

These formulas have been derived assuming that the velocity and buoyancy 
orce are constant across the jet and zero outside of it and assuming that the 
et had no initial momentum. 

& 
aaa 
9 «Turbulent Gravitational Convection from Maintained and Instantaneous Sources,” 


B.R. Morton, G. Taylor and J. S. Turner, Proceedings, Royal Soc. of London, Series 


., March 6, 1956, 11, pp. 1-23. 
10 «Gravitational Convection from a Boundary Source,” by H. Rouse, C.S. Yih and H. 


v7. Humphreys, Sartoyck ur Tellus, Vol. 4, No. 3, August, 1952, pp. 201-210. 


6 June, 1960 HY | 


The formulas are only valid for a jet satisfying d= 0. For a jet with finite 
diameter the factor x of Eqs. 14 and 15 must bereplaced by the factor (x + 2d) 

The constant @ must be found by experiment. Experiments made with wate: 
jets indicated that a = 0.093. 

It is possible to write the equivalent of Eqs. 14 and 15 for a jet with a finit 
diameter ina dimensionless form by use of the relationships in Eqs. 4 an 
a@ = 0.093. This gives 


1/3 
am = 3.05 F2/3( 3 Re Cee (16 
u d 
Oo 
and 
c ARES: 
fm _ gq p2/3 (2 +2) 6 ee al (17 
Co d 


The experiments of Morton, et al.,9 showed that the distribution perpendicu 
lar to the axis of the jet is given by the Gaussian curves: 


and 
C= Cin 


Rouse, et al.,10 also derived the 5/3 and 1/3 power of Eqs. 16 and 17 theo 
retically. Their semi-empirical formulas, based on experiments done with ai 
jets reduce in dimensionless form to: 


Um _ -2/3 [x\- 1/3 
ae os) 4.35 F (3) eo 0 we. ape ee eee (2( 
= 5/3 
m BY Blox 
co 7 935 F (3) One ee (2 
| - 96(y/x)2 
We UO a ot eA Se aon eee (28 
and 
es - T(y/x)? 
CrP CCE eS > eS. iy as 2 Ra es Peau 


The author’s experimental results agree better with Eqs. 16 and 17, due t 
Morton, et al., than with Eqs. 20 and 21, due to Rouse, et al., since the firs 
set of equations takes a small initial momentum in account. Therefore th 
author has preferred to apply Eqs. 16 and 17 in the following section. 7 

The Intermediate Case, Finite Values of F.—The intermediate case has bee 
studied theoretically by Schmidt.3 The equations that he derived suggest the 
the intermediate case may be described approximately by the equations for th 
non-buoyancy case (Eqs. 6, 7, 8, 9, 10, 11 and the previously selected value 
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fk, uw, Cy, Co, C3 and C4) from the nozzle to a certain distance from the 
nozzle, and by the equations for the buoyancy case (Eqs. 16 and 17) at greater 
distances from the nozzle. The transition occurs at the point x = xX, where 
the curves that describe the non-buoyancy case and the buoyancy case inter- 
sect each other (Fig. 1). 

For sufficiently large values of F (F > about 10) the ordinates x; of the two 
sets of points of intersection satisfy 


2Co 1 /%\? 
ST a 0 eee Alic36 DESEO (24) 


with Cy, 4. and k according to the previously selected values. 


For sufficiently large values of x; and hence for sufficiently large values of 
F 


n d2 
Ane fo) se 


2 Cy 4 ( zy - ef eliss(-0m 


since for sufficiently large values of x; the contribution of the zone of flow es- 
tablishment to the integral of Eq. 25 is relatively small and then it does not 


make much difference if we apply Eqs. 10 and 11 also for a <6 or 7 to evalu- 


ate the integral. 

Eqs. 24 and 25 show, for sufficiently large values of F at least, that the in- 
fluence of the initial momentum is larger than the influence of the buoyancy 
force for x< x;, while the contrary applies for x > x. This is compatible with 
the description of the intermediate case suggested by the equations of Schmidt 
of which Fig. 1 is a graphical representation. In the following section experi- 
ments will be described by which the curves of Fig. 1 are confirmed. 

Fig. 1 only applies for the case that the initial momentum and the buoyancy 
‘orce have the same direction. For the case that the initial momentum and the 
uoyancy force have opposite direction, the distance x; is an approximation of 
he distance from the nozzle, where the velocity along the axis will be reduced 
O zero. 


EXPERIMENTS 


The experiments to check the families of curves of Fig. 1 were made ina 
rlass-wall tank, with a height of 3 ft and a cross section of 1 ft by 1 ft. The 
et issued into the tank through pipes through the tank bottom with internal di- 
meters of 0.6 in., 0.8 in. and1.0 in. The end of the vertical pipes was about 
| in. above the bottom of the tank. The upper edge of the tank was constructed 
1s a continuous weir. The discharge of the jet was determined by a volumetric 
neasurement of the discharge over the weir. 

The density of the fluid of the jet was about 1,000 kg per cu m. The density 
f the surrounding fluid varied from 1,020 kg per cu m to 1,050 kg per cu m. 
the surrounding fluid was a solution of crude table salt and water. 

Measurements of the Density.—The density at the axis of the jet was mea- 
ured to check the curves for the soncentzation: Samples of the fluidof the jet 
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FIG, 2; —EXPERIMEN TAL SETUP 
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were taken at different points in the jet through pipes with an internal diameter 
of 3/16 in. The pipes were mounted to a frame that could be moved through 
the fluid in such a way that the inlet opening described a vertical line (Fig. 2). 
The discharge through the pipes was such that a periodof 25 sec was required 
to collect a sample of 100 cc. The density of the samples was measured with 
a Westphal balance. The balance was accurate to about 0.1 kg per cu m. 


TABLE 2.—RESULTS OF MODEL RUNS 


p Theory 
Run he in x sae wie F R <m <m events 
inches d =n Lae cS Co in % 
theory of Co 
21 0. 16.4 2. 0.235 1 ly 
23 0. 19.2 2. 0.10 +2 
15 0. 16.8 2s 0.175 -1/2 
16 0. 21.0 2. 0.08 ctaely2) 
id 0. 18.1 2. 0.085 0 
22 0. 21.2 3. 0.15 Ad 
8 0. 22.0 3. 0.215 + 1-1/2 
3 0. 20.2 3. 0.22 0 
6 0. 27.2 3. 0.33 + 1-1/2 
10 0. 30.8 3. 0.23 +3 
3 0. 4, +1 
5 0. 5. om) 
12 0. 5 0 
a 0. 5. 12 
2 0. 5 -1/2 
7 0. 10 Wp -6 
46 0. 26 7 0 
43 0. 23 8 cae ath 
44 0. 26 9 + 1-1/2 
36 0. 4 1 +3 
38 0 6 i 0 
35 0 4 2 -9 
42 0. 1 2 +4 
47 1 6 1 ol 
51 1 0. 1. 72 
48 1 6 i 
52 a 0 AL 
54 1 2a ne 
49 a 6. 2 
53 it 0. 2 
50 1 6. 2 


_ Thetime history of the density of the fluidof the jet was measuredat differ - 
wnt points in the jet. A typical example is shown in Fig. 3. Fig. 3(a) shows 
he time history on ordinary scale and Fig. 3(b) shows the same data on asemi- 
ogarithmic scale. The first steep part of the logarithmic plot (and the corre- 
ponding part of the other plot) describes the establishment of equilibrium. 


é 
. 
4 


: 


12 June, 1960 HY 6 


The second less steep part describes how the equilibrium is disturbed by the 
fact that the surrounding fluid is gradually mixed with the fluid of the jet. The 
small dimensions of the tank made it impossible to obtain a real equilibrium. 
Therefore, the density corresponding to the point of intersection of the two 
straight lines of the logarithmic plot has been considered to be the equilibrium 
density. Combined with any error made in the density measurement itself, 


0.8 
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FIG, 4.—COMPARISON BETWEEN THE EXPERIMENTAL 
RESULTS AND THE THEORETICAL DATA 


this method will cause an error of about 0.5 kg per cu m in the determinatior 
of the density difference between the fluid of the jet and the surrounding fluid 


- The dilutionis expressed in percentage. It is defined as Sn 2 nc Weta x 100 


: ; Obs Ps =P 
For the experiments, p, - Po Varied between 20 kg per cu m MAL 5 per ci 


m. Hence the possible absolute error in the determination of the dilution is 


O20 yee 
about 30 to 50 2%. 
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The determining circumstances and the results of the model runs regarding 
the density are shown in Table 2. For most of the runs the Reynolds number 
is sufficiently high to assure turbulent flow conditions. 

The experimental results are plotted in Fig. 1. A comparison between the 
experimental results and the theoretical data, based on the theory of Morton, 
et al.,9 is shown in Figs. 4(a) and (b). 

The experimental results agree well with theory. The approach, suggested 
by the equations of Schmidt, to combine the results of the non-buoyancy case 
with those of the buoyancy case yields a good description for the intermediate 
case. That Eqs. 16 and 17, due to Morton, et al.,9 take a small initial momen- 
tum in account is the reason for the goodagreement obtained for small values 
or F. 

Measurements of the Velocity.—A pitot tube has been used in an attempt to 
measure the velocity-distribution in the jet. The velocity readings were, how- 
ever, insufficiently accurate because of the following reasons: (a) equilibrium 
conditions were not reached in the tank during the model runs; (b) the velocity 
in the jet is a function of the distance from the nozzle, while a pitot tube re- 
quires a constant velocity between the static- and dynamic-opening; and (c) the 
small inaccuracies made in the determination of the velocity-head were im- 
portant because of the magnitude of the velocities involved (20-30 cm per sec). 

Since no other instrument (small propellor, for example) was available, no 
further attempt to measure the velocities has been made. 


CONCLUSIONS 


The family of curves shown in Fig. 1(a) describes the dilution due to a sub- 
merged circular turbulent jet, issuing vertically upward in a fluid with a higher 
specific gravity than the jet. The curves are based on theory, which is ade- 
quately confirmed by model experiments. The density-range covered by the 
experiments is such that the curves may be used when sewage disposal in the 
marine environment is considered. 

The family of curves, shown in Fig. 1(b), describes the velocity -distribution 
in the jet. The curves are based on theory. No attempt has been made to con- 
firm this set of curves experimentally. 
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TRANSLATIONS OF FOREIGN LITERATURE ON HYDRAULICS 


Second Progress Report of the Task Force on List of Translations of the 
Committee on Hydromechanics of the Hydraulics Division 


INTRODUCTION 


In an attempt to bringto hydraulic engineers an up-to-date listing of trans- 
lations of foreign literature on hydraulics, it is the intent of the task force pre- 
paring this report to issue, fromtime totime, progress reports containing suit- 
able items. Interested readers are urged to submit discussions to (a) add to 
this list, (b) offer suggestions for improvement, and (c) otherwise assist the 
task force infulfilling its aims. The present listis the second issued since the 
publication, in 1957, of ASCE Manual 35 entitled “A Listof Translations of For- 
eign Literature on Hydraulics.” Thefirst list appearedin the August 1959 Jour- 
nal of the Hydraulics Division. This list is to be considered as an addendum 
to Manual 35 and will be includedin a revision of that Manual when the time is 


deemed appropriate. 


ABRAMOVICH, G. N. “The Turbulent Jetin a Moving Fluid.” Izv. Akad. Nauk, 
Otd. Tekh, Naul. 1957, No. 6, pp. 93-101, M. 465. (On file at the Dept. of 
Scientific and Industrial Research, Lending Library Unit, London, England). 


AL’ TSHULER, L. V., BAKANOVA, A. A., AND TRUNIN, R. F., “Phase Trans- 
formations During Compression of Water by Strong Shock Waves.” Akade- 
miya Nauk SSSR. Doklady, Vol. 121, No. 1, pp. 67-69, 1958. Tr. No. 59-15750. 
(On file at LC or SLA). 

ANONYMOUS. “Symmetrical Cavitation Flow Arounda Flat Plate, Situated Be- 
tween Parallel Walls.” Akademiya Nauk SSSR. Otdeleniye Tekhnischeskikh 
Nauk, Izvestiya, No. 4, 1946. Tr. No. 59-15095. (On file at LC or SLA). 
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_ Note.—Discussion open until November 1, 1960, Separate discussions should be sub- 
mitted for the individual papers in this symposium. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 86, No. HY 6, June, 1960. 
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ANONYMOUS. “Instrument for Measuring Cavitation During the Operation of 
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Izobretatel’stvo V, SSSR, Vol. 2, No. 10, p. 40, 1957, DSIR 33866/CT. (On 
file at Dept. of Scientific and Industrial Research, London). 


ARKHIPOV, V.N., “The Formation of Streaming Fluctuations Behind a Solid 
Obstacle.” Dok Ak Nauk SSSR, Vol. CXXII, No. 4, pp. 620-624, 1958. Tr. 
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FLOW LOSSES IN THE LOWER GILA RIVER 


By Lawrence F. Pratt,! F. ASCE 


SYNOPSIS 


This paper describes a method for estimating losses from infrequent flows 
in the lower reaches of the Gila River, in Arizona. The same procedure could 
probably be applied to other intermittent streams if sufficient data are avail- 
able. 


INTRODUCTION 


The usual procedure for studying or analyzing flood flows in a stream 
wherein water flows continuously is some kind of mass-curve technique. The 
literature concerned with this problem is voluminous. However, such proce- 
dure has little application to intermittent streams because the available data 
are too discontinuous to provide a suitable universe for common statistical 
methods. 

Because of the discontinuous nature of the gaging records, the only ration- 
al way to study flood flows of an intermittent stream is to take them one at a 
time, not enmasse. At a gaging station on an intermittent stream, there will 
be periods of weeks, months, even years, when no flow will be recorded. Be- 
cause of the many zeros in the record any sortof average is meaningless. Any 
given zero flow may reflect a condition such that the last surface water may 
have disappeared only a few feet, or many miles, above the gage. If there is 
another gaging station upstream, it is obvious that the upper gage will record 
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many more daily flows than the lower gage, provided that local inflow between 
the gages is negligible or absent. 

It is possible to express the percentage of loss between an upstream and a 
downstream station interms of the quantity of flow at the upstream station and 
duration of flow at the downstream station. The procedure is to cumulate con- 
current flow quantities at both stations, making proper allowance for time of 
travel between the two. This allowance can most conveniently be made by 
matching the days of peak flows, but is not constant. When the streambed has 
been dry for a long time, much water will be required to fill surface depres- 
sions and to initiate infiltration. These requirements take time as wellas wa- 
ter. When flow has been established, the time of travel may be reduced to as 
little as a third of the time of travel for a dry streambed. 

As flow continues, the percentage of water lost betweenthe upper and lower 
stations changes, rapidly for the first few days and then more slowly, ap- 


fe) 
PHOENIX 


SAN 
Ce, o CARLOS 


° 
FLORENCE 


FIG, 1.—LOCATION MAP 


proaching a limit. If flow continues long enough to indicate the value of the 
limit, it may be inferred that this value is primarily caused by evaporation, 
since equilibrium will have been approached between surface flow and infil- 
tration. Higher percentages of loss represent the prior mentioned require- 
ments for filling depressions and initiating infiltration. | 


GILA RIVER BASIN 


An important intermittent stream of southwestern United States is the low- 
er portion of the Gila River, in Arizona (see Fig. 1). The reaches considered 
herein are from Dome to Gillespie Dam, since 1928, and from Dome to Mc- 
Dowell and San Carlos, for 1904 and 1905. | 

The Gila River Basin covers an area of 58,180 sq miles, of which about 
88% is in Arizona, 10% is in New Mexico, and 2% is in Sonora. The upper 
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reaches of the Gila River and its principal tributaries are in mountainous 
terrain within which are several peaks of over 10,000 ft elevation. It is from 
this terrain that nearly all runoff is derived. At about 250 to 300 miles above 
the mouth of the Gila River, several storage dams have been constructed to 
impound this runoff. The larger dams are: Coolidge on the Gila; Roosevelt, 
Horse Mesa, Mormon Flat, and Stewart Mountain on the Salt; Horseshoe and 
Bartlett on the Verde; and Carl Pleasant on the Ague Fria. In 1952, the usable 
storage capacity of reservoirs on the Gila River and its tributaries was 3,446, 
000 acre-ft2. This storage capacity is not quite twice the average annual in- 
flow to the reservoirs. 

Downstream from the reservoir system is the Salt River Valley contain- 
ing some 800,000 acres of irrigated land and several cities and towns, the 
largest of which is Phoenix. In nearly every year, all available surface wa- 
ter released from storage or running in unregulated tributaries is wholly con- 
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FIG, 2.—FLOW AT GAGE FROM JULY 30, TO NOVEMBER 26, 1904 


sumed by irrigation and domestic uses. In addition, a large quantity of ground 
water is pumped and consumed. The result is that very little water appears 
below Gillespie Dam, the outlet of the valley, and only rarely does any water 
reach the mouth of the Gila River or the gaging station at Dome, 12 miles 
above the mouth. The river distance from Gillespie Dam to Dome is about 
150 miles. Many of the flows registered at the Dome gage originate from in- 
frequent rains falling west of Gillespie Dam or an return flows from irrigated 
areas between the two gages. 

_ Early Records?.—For a brief period during the years 1904, and 1905, gag- 
ing stations were operated concurrently on both the Salt and Verde Rivers just 
above their junction at McDowell, and on the Gila River at San Carlos (near 
Hi : 


5) 2 «Ground Water in the Gila River Basin and Adjacent Areas—A Summary,” An open 
e report of the USGS, 1952. 
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present Coolidge Dam) andat Dome, also called Gila City. During this period, 
there were-no storage reservoirs in the basin. 

The Gila River was dry at Gila City (Dome) from January 1, 1904, (and for 
an unknown time earlier) to July 29, 1904. Flow at the gage started July 30, 
and continued until November 26, after which date the river was again dry. 
Only gage heights are published for Dome, but a rating curve was prepared by 
the author, using data for other periods, to estimate the average daily flows. 
Average daily flows are published for the three upstream stations. 

The apparent time of travel of flood-peaks, after flow was established at 
Dome, was seven days from McDowell to Dome, and five or six days from San 
Carlos to Dome. While the flow passing Dome persisted for 119 days, at no 
time did it become large. The maximum daily flow at Dome was 4,560 day- 
sec-ft and the total flow was about 220,000 acre-ft. The combined flows pass- 
ing the three upstream stations was about 508,000 acre-feet. The difference, 
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FIG, 3.—FLOW FROM JANUARY 16 TO JULY 19, 1905 


or apparent loss, was 288,000 acre-ft. This difference ignors both unmeasur- 
ed inflows and consumptive uses, neither of which can be estimated with any 
assurance for this period. 

Flows at McDowell and at San Carlos were combined at six days earlier 
than Dome. The combined flows were cumulated day by day and the same cu- 
mulation was made for the flows at Dome. The cumulative differences were 
considered to represent the cumulative losses and were tabulated both in day- 
sec-ft and as percentages of the combinedupstream flows. Percentage of loss 
to the end of the flow period was 56.8%. Cumulative loss percentages betweer 
55% and 57% were consistent for more than the last half of the period (seeé 
Fig. 2). It is concluded that this order of loss percentage fairly represent: 
the conditions existing in 1904, with a relatively low but persistent flow. : 

The Gila River was dry at the Dome gage from November 26, 1904, unti 
January 16, 1905. Flow began on that date and continued until July 19, 1905 
except for February 2, 3, and4. Very large flood peaks occurred in this si: 
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months period, reaching maxima of 82,000 day-sec-ft on February 8, 95,000 
day-sec-ft on March 20, and 64,000 day-sec-ft on April 14. Except for a re- 
latively small peak of 13,000 day-sec-ft on April 28, flows gradually decreas- 
ed from April 14, to July 19, when the river was again dry at Dome. Because 
of the numerous large flood peaks, the calculated cumulative percentage of 
loss is not as constant as was that calculated for the 1904 period (see Fig. 3). 
For the last half of the period the loss was fairly steady at around 40%, but 
during the first half the loss varied from 21% to 43%. The measured inflow 
for the period was about 5 million acre-ft and the outflow at Dome was about 
3 million acre-ft. Again, consumptive uses and unmeasured inflows were 
ignored. 

Later Records¢.—From 1907, until a recording gage was established in 
May, 1929, there are no records of discharge at Dome. A recording gage had 
previously been established at Gillespie Dam in July, 1924. From May, 1929, 
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DURATION OF FLOW - DAYS 


FIG, 4.—GILA RIVER, 1929 TO 1955 


through September, 1955, there were a total of 312 months. In only 29 of those 
months were there recorded flows at Dome in excess of 1,000 acre-fta month. 
For the same period there were 90 months when the flow past Gillespie Dam 
xxceeded 1,000 acre-ft. Therefore, two-thirds of the months when there was 
appreciable flow past Gillespie Dam there was no flow at Dome. Under these 
‘ircumstances, a correlation of monthly flows at the two stations is quite 
meaningless and a correlation of annual flows is impossible. 

During the 27 yr (1929-55) of records studied on the Gila River, there were 
mly eleven occasions when flows passing Gillespie Dam reached the Dome 
sage for sustained periods of seven days or more. The longest periods of sus- 
ained flow were in 1941, for 114 days, and in 1932, for 66 days. Two flows, 
oth in 1931, continued for 22 days and 25 days. Seven other flows were sus- 
‘ained for 7 days to 12 days. 


a 


_ 4 USGS annual water supply papers for the Colorado River Basin. 
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It was found that the percentage of the accumulated flow passing Gillespie 
Dam, but lost before reaching Dome, was related to the duration of flow at 
Dome. It was also found that the time of travel was affected by the antecedent 
condition of the streambed. If the bed had been dry for several months, the 
travel time was aS much as eleven days. After the channel had been wetted, 
the time was reduced to three or four days. The time of travel for each per- 
iod or sub-period was determined by matching the appearances of flood peaks 
at the two gaging stations. The differences between the cumulative flows at 
Gillespie Dam and the cumulative flows at Dome were found and expressed as 
percentages of the Gillespie flow: These differences were considered to be 
the cumulative losses in the reach, ignoring any inflow caused by local pre- 
cipitation. The determination of any such local inflow would involve the rough- 
est kind of estimates, as there are no records of flow and an insufficient num- 
ber and inadequate distribution of rainfall records to make useful rainfall- 
runoff correlations possible. Consumptive uses were small and were ignored. 
The results were plotted as Fig. 4. 

Fig. 4 shows that the percentage of loss drops quickly for the firstfew days 
and then flattens out, usually approaching a limit, sometimes gradually anc 
sometimes quite abruptly. For the two flows of more than 30 days duration 
the curves are of hyperbolic form, which suggested the calculation of an em- 
pirical equationto represent the average of these two occurrences. The equa- 
tion is of the form: 


in which y is cumulative percentage loss, ais the first day percentage loss 
x is days of flow at the lower gage, and b is an empirical exponent. For the 
reach from Gillespie Dam to Dome the equation is: 


Se ered a ee eee (2 


CONC LUSION 


It seems probable that an equation of the same form could be derived fo: 
any other intermittent stream if sufficient data were available. The constant; 
would undoubtedly be different. 

For those streams wherein flows are more frequent and of longer averag 
duration than for the lower section of the Gila River, a relationship linkin 
duration of flows to channel losses, or to flows at the lower station, woul 
have more significance. If the economic value of the water could be establish 
ed it would then be possible to determine to what extent salvage should b 
undertaken, by storage or channel control, to increase the duration of flow 
and thereby decrease the cumulative percentage of loss. 
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CAPACITANCE METHOD OF MEASURING WATER FILM THICKNESS2 


By R. H. Black! 


SYNOPSIS 


A water film thickness gage was developed for making rapid measurements 
of the depth of a water sheet flowing ona plane surface, The gauge uses the 
electrical capacitance between a fixed metal plate and the water surface to 
sense the air gap and, by use of differences, measure the water depth. 

Rapid readings of the depth of a disturbed water sheet flowing down a plane 
were made, The dynamic range of the instrument was 0,130 in. with an accu- 
racy of +0.003 in. Transverse roller waves 3 in. long were resolved with this 
instrument; through the use of a special “probe,” waves 1/2 in, long were re- 
solved. 


INTRODUCTION 


The U. S. Naval Radiological Defense Laboratory is currently studying the 
transport of radioactive particulate contaminants by water sheets, Stream depth 
is one of the most important factors in the equations describing sheet-flow of 
water on an incline plane. 2 This depth can be measured precisely and accu- 


EL Re a a ee a oan ee 
Note.—Discussion open until November 1, 1960. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 86, No. HY 6, June, 1960. 

_ 2The opinions expressed by the author do not necessarily reflect those of the Navy 
Dept. : 

_ 1Chem. Engr., Aerojet-General Nucleonics San Ramon, Calif, Formerly with the U. 
§. Radiological Defense Lab., San Francisco, California. 

_ 2Transport of Contaminant by a Water Film,” by R.H. Heiskell, R. H. Black, and 
H. L,. Burge, U. S. Naval Radiological Defense Lab. Review and Lecture, USNRDL-RL- 
36, February, 1958. j 
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rately to +0.001 in. by the use of a micrometer point gauge when stream flow 
is uniform and the surface is glassy smooth. However, when water flows down 
an incline plane with a slope greater than 0.01, pronounced transverse roller 
waves are generated a short distancefrom the entrance to the plane, and these 
waves introduce serious problems in the use of a point gauge for measuring 
depths, Wave height is not a constant for a given flow condition and trough 
measurements are questionable, as the “end point” for determing contact is 
subjective. Furthermore, no information concerning wave-shape results from 
the use of this techique. 

A review of the literature revealed two other methods of determining water 
sheet thickness. These were an optical method,’ which requires the water to 
be flowing on the outside of a cylinder, and a capacitance method,‘ utilizing a 
frequency-modulated carrier, The first method was impossible due to the un- 
suitable geometry of the apparatus required tobe usedand the second was con- 
sidered likely to be subject to electronic drift. 

As future progress in the experimental study of sheet flow depends partly 
on the development of an instrument that can be used for making depth meas- 
urements during non-uniform flow, the objective of this project was to develop 
an instrument capable of making depth measurement of a sheet of water flow- 
ing over a plane surface, and to obtain adequate response to demonstrate wave 
shape. 


APPROACH 


There were several potentially useful means for determining the thickness 
of a water sheet, grouped under two basic approaches. The first was to utilize 
absorptive properties relating to the thickness of the water sheet by (a) de- 
tecting the absorption of beta radiation passing through the water, or (b) de- 
tecting the absorption of light passing through dyed water. The second basic 
approach was to measure the distance to the air-water interface and obtain the 
thickness by ( a) subtracting the distance to the stream bed by detecting con- 
tact with the surface electrically or optically, (b) detecting the beta backscatte1 
of a radioactive source fixed above the surface, or (c) detecting the electrica! 
capacitance between a metal plate fixed above water and the (electrically ground- 
ed) water surface. 

A brief consideration of these possible means from the viewpoint of utiliz- 
ing readily available equipment revealed that most of them had shortcomings 
Methods dependent on radioactive sources would require relatively expensive 
equipment and heavily shielded radiation sources, Electrical and optical-contac 
detection equipment was notavailable, Light-absorption methods, utilizing dye 
water, would necessitate a recirculating system anda translucent or trans: 
parent stream bed, and would complicate or delay other studies requiring 
transparent water sheet, However, a capacitance transducer was commerciall: 
available utilizing a patented mechanic-electric transducer which had a spec. 
ified sensitivity of 2 v output per ufd change of capacitance and less tha 
0.1% drift per.° The detection of electrical capacitance between a fixed ref. 

3 “Heat Transfer by Condensing Vapors,” by C. G. Kirkbride, Transactions, Amer 
Inst. of Chem, Engrs., 30:170, 1934, 

4 “Characteristics of Flow in Falling Liquid Films,” by A. E. Dukler and O. P. Ber- 
gelin, Chem. Engr. Progress, Vol. 48, No. 11, 1952, p. 557. 
ae “Mechanical - Electric Transducer,” by Kurt S, Lion, Rev. Sci. Instruments, 27:222 
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rence plate and the water surface would avoid the disadvantages of contact 
ith the water or dyeing of the water. 


APPARATUS 


The commercially available capacitance transducer is essentially a capac- 
ance comparator and gives adirect current (D. C.) signal that is proportional 
the percentage difference between the values of two capacitors witha common 
ound.© The transducer is excited by and permanently connected to a 250 kc 
rier generator. One of the capacitors, C+, is a plate, or the “probe,” fixed 


ove and parallel to the electrically grounded water. The other, C“, is an ad- 
stable reference capacitor. The relative positions of these and the other elec- 


ical units are shown in Fig. 1. They are described in the Appendix. The R. 
filter and D. C. power supply were constructed according to specifications 
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FIG, 1.—BLOCK DIAGRAM OF INSTRUMENT 


the manufacturer of the capacitance transducer. The D. C. amplifier and 
adout can be either anoscilloscope with acameraor a recording oscillograph. 
The probe area was arbitrarily chosen as 1.0 sqin. (6.42 sqcm). Capacitor Cy 
designed to vary +0.3ufd with water-sheet thickness changing by 0.130 in. 
330 cm), and to have a mean capacitance, including fixed stray capacitance, 
3 to 5yfd. Capacitor Cy is adjustable from 3 to 12pfd. The capacitance 
tween the probe and the water surface is considered to be the same as that 
‘ween two parallel plates as a first approximation in setting the probe para- 


ters; 


= 0.0885 © eres ¢-etent 1) 


C, : - andtecbgenceria le 


which C1 is the capacitance (ufd), s denotes probe area (cm), and d is 
> distance from probe to water (cm) 
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With the values given previously, two equations can be established: 


6.42 


0.0885 x Pd? ets ee eae hae sect aoeee ( 


iH} 


Cy 


x 6.42 ( 


0.0885 (d + 0.330) cospsrpveecuun se Que 


ii} 


Cy + 0.3 


Solving these equations gives 0.6pfd for Cy and 0.38 in. (0.97cm) for d. 
The probe (which is removeable) and a shield, that encloses and fixes t 
stray capacitance of the transducer, reference capacitor, and probe lead wir 


DIFFERENTIAL 
C.AMPLIFIER 


FIG, 2.—APPARATUS 


are rigidly fastened to the probe mount. This assembly has been offered 
cently as a unit by the transducer manufacturer. The probe mount is held 
vise which rides on a vertical calibrated micrometer screw. With this 
crometer screw, the probe can be raised or lowered by a known distance f 
the stream bed. The instrument is shown in Fig. 2. 
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CALIBRATION AND PERFORMANCE 


A reasonable compromise between adequate sensitivity and an acceptable 
lynamic range was found with a probe 1-1/4 in, wide by 3/4 in. long, designated 
is the “standard” probe. With the standard probe set 0.400 in, above the stream 
ed, and the differential D. C. amplifier gain set to give full-scale response to 
.1-v signal, the dynamic depth range is 0.030 to 0.160 in. 

The maximum sensitivity of the instrument for measuring change of depth 
S approximately 0.0003 in. This sensitivity can be used for measuring water 
urface undulations; however, accuracy is not great enough to use the full sen- 
itivity effectively in determining water film thickness. The sensitivity of the 
nstrument utilizing the 0.130 in. dynamic range is +0.001 in, 

The instrument was calibrated in place with the water flowing. A 4 in,-by- 
-in. metal plate of known thickness (and approximately the thickness of the 
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FIG, 3.—TYPICAL RECORDINGS ILLUSTRATING 
CALIBRATION AND OPERATION STEPS 


ater sheet) was placed under the probe to displace the water and give a mid- 
cale reference reading. The probe was lowered by advancing the micrometer 
crew until the output meter barely read off-scale; then the probe was raised 
y 0.010 in. increments (or other convenient mena until the opposite scale 
mit was reached, and incremental readings were recorded to provide a scale. 

ollowing this, the probe was re-set to the mid-scale reading, and the reading 
r the plate was recorded to provide a reference, thus, completing the cali- 
ration, The reference plate was removed, and the instrument was then meas- 
cing the water sheet thickness. Fig. 3 shows a specimen ofthe calibration 
ta and depth data. 

The accuracy of the instrument was determined by comparing its readings 
ith water depth measurements of smoothly flowing water taken with a mi- 
ometer point gauge. Agreement was + 0,003 in. 

x 


: 
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The resolution of the instrument to longitudinal waves using the standar 
probe was measured by replacing the water sheet with a flat metal plate, Thi 
plate had several widths of slots cut to adepth of 0.030 in. to represent “waves. 
This plate was moved under the probe, in place of the water sheet, and the slc 
depth, as indicated by the instrument, was recorded. The results of this tes 
are given in Fig. 4. A 0.5 in. slotcan be resolved with a 50%: response to ampli 
tude, and a 1.5 in. slot isthe shortest that can be resolved withnoloss in ampli 
tude. This resolution is assumed to be the same as the instrument resolutio 
to longitudinal water waves. 


0.100 


0.030 “Se 


0.010 


0,003 


INDICATED DEPTH (IN) 


AIR GAP BETWEEN 
PROBE AND PLATE 


0.015 IN 
0,100 IN. 
0.200 IN. 
0.300 IN. 


SLOT WIDTH (IN.) 


FIG, 4.—RESPONSE OF INSTRUMENT TO 0.030-in, 
DEEP SLOTS IN A METAL PLATE 


A special probe was constructed to resolve shorter waves. This small pro 
is essentially a short 0.080 in. diameter rod fixed perpendicularly tothe strez 
bed and shielded along its length with a 0.25 in. diameter tube. The results 
resolution tests for this smaller probe are also presented in Fig. 4. A 0.1: 
slot can be resolved with a 50% response to amplitude. This probe has li 
sensitivity and very nonlinear response to water-sheet thickness, and, therefo 
was not generally used for thickness measurements. However, it can be ca. 
brated in the same manner as the standard probe. . ‘ 


CONCLUSIONS 
The accuracy of the instrument, as determined by comparison with the n 
crometer point gauge, is limited by the conductivity of the water. Although 1 
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mpedance of the transducer is high relative to that of the water, slight fluc- 
uations in the conductivity of the water due to depth and ion content affect the 
esistance and, therefore, the apparent “location” of the grounded side of Cy. 
imple tests indicated that accuracy canbe improved to +0.001 in, by increas- 
ng salinity to approximately 0.01 Normal NaCl. = 

The resolution of the instrument to waves is, in part, limited by the dimen- 
ions of the probe. The shortest logitudinal wave that can be resolved is equal 
0 the length of the probe. 

When the small probe is used, the resolution of waves is limited by diver- 
ence of the electrostatic lines of force from the smallarea of the probeto the 
ater surface. If the electrostatic lines could be made parallel by applying 
he proper voltage to an electrostatic ring around the probe, then, the resolu- 
ion of the instrument would be limited only by the smallest probe which could 
upply a readable signal. 

The instrument sensitivity to changes in water depth is limited by the elec- 
rical noise level (which is approximately 2 millivolts), the maximum probe 
rea as determined by wave resolution requirements, and the minimum air gap 
etween the probe and the water surface, as determined by the anticipated fluc- 
lation in depth. 

Instrument stability is excellent; drift was within +0.001 in. over a 4-hour 
eriodaftera 15 min. initialwarm-up. The transducer life appears reasonable; 
le unit required repair after one year of modest use. 

Instrument sensitivity, useful range, and wave definition are interrelated, 
s indicated above, must be balanced against each other for a particular appli- 
ation. By manipulating the value of Cg, the area and shape of the probe, the 
istance of the probe from the mean water surface, and the type of readout or 
ecording equipment, this depth-measuring instrument can be applied to a va- 
iety of hydraulics studies involving afreeliquid surface. The rapid indication 
f water surface and ease of reading, even with a disturbed surface, can save 
uch time. Theavailability of electrical signal output suggests the possibilities 
‘centralized simultaneous recording of multiple stations andautomatic control. 


APPENDIX 


Equipment used for this report: 
apacitance Transducer; Delta Unit Model 901-1, Decker aviation Corp., Phil- 
lelphia, Pa., 2V/ufd sensitivity at 10 ufd initial capacity. 


egulated Power Supply; NRDL design (vacuum tube regulation) to comply with 
ecifications furnished by manufacturer of capacitance transducer. 

ilter; NRDL-fabricated per diagram supplied in capacitance transducer in- 
ruction manual, 


irect-writing Oscillograph; Sanborn, Model 128, 200 mv/cm sensitivity, dif- 
rential D. C. input. 


scilloscope; Dual Beam, Type 322A, A. B. DuMont Lab., Clifton, N. J. 200 
v/em deflection sensitivity, differential D. C. input (one beam required). 


scillograph-record Camera; Type 297, A. B. DuMont Lab., Clifton, N. J. 
mstant Voltage Transformer; Solavolt No, 50106, 500 Watt. 
y 
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SEDIMENTATION ASPECTS IN DIVERSION AT OLD RIVER 


By Fred B. Toffaletil 


SYNOPSIS 


In recent years the Old and Atchafalaya Rivers have diverted a progres- 
vely increasing percentage of Mississippi River flow. The damaging conse- 
ences posed by this continuing action necessitated controlled diversion. This 
per presents the sedimentation aspects of channel evolution leading to ne- 
ssity of control. It is shown that the distributary system is in a state of im- 
lance in sediment transport capacity, the nature of which led tothe conclusion 
at the flow capacity of the Atchafalaya River canbe expected to increase annu- 
ly regardless of the type of water year. This conclusion is substantiated by 
ld observations at the Simmesport discharge range. 


INTRODUCTION 


The Atchafalaya River, which is connected with the Mississippi River by Old 
ver, is the farthest upstream distributary of the Mississippi River. Fig. 1 
a plan map showing the Mississippi River below latitude Natchez, the lower 
J of the Black and Red Rivers, the Atchafalaya River, and connecting Old 
ver. From the head of Old River it is approximately 140 miles to the Gulf 
Mexico by way of the Atchafalaya River. By way of the Mississippi River 
3t Baton Rouge and New Orleans it is 320 miles, or more than twice the dis- 
.ce via the Atchafalaya River. Old Riveris approximately 6.5 miles in length. 


Note. —Discussion open until November 1, 1960. To extend the closing date one month, 
ritten request must be filed with the Executive Secretary, ASCE. This paper is part 
copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
a Civil Engineers, Vol. 86, No. HY 6, June, 1960. 
Chf., Hydr. Model and Sedimentation ‘Sect., Hydr. Branch, U. S. Army Engr. Div., 
ver Missippi Valley, Vicksburg, Miss. 
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The lower 60 miles of the Red River channel and the lower portion of Blacl 
River, whose principal flow is from the Quachita River, are in a backwate: 
basin containing numerous shallow lakes and low-level swamp areas. At thi 
peak of Mississippi River floods, and coincident peak flows in the Atchafalay: 
River, the water-surface slopes through this backwater area are very flat 
Consequently, the Red and Ouachita Rivers system can contribute only ver 
fine sediments to the Atchafalaya River at high stages. 

In recent years the Old and Atchafalaya Rivers have diverted progressivel 
increasing percentages of Mississippi River flow. Studies of this developmen 
led to the conclusion that this diversion from the Mississippi River will reac 
a “critical” stage (approximately 40% of the major river’s flow) by 1975. Th 
consequences of such an eventuality necessitated a controlled diversion. Con 


FIG, 1.—RIVER SYSTEM BELOW LATITUDE NATCHEZ 


struction of control structures was initiated in 1955 and is now essential 
complete. All work, including new channels, bank protection, a lock, and tt 
closure of Old River, is scheduled for completion by the end of 1963. Detai. 
of this project which include river history, the general problem, hydraulic re 
quirements of control works, type structures and foundation design, have bee 
presented previously. 

Previous discussions of this problem area have dealt mostly with flow a1 
channel dimensional changes. This paper presents the sedimentation aspec 
of channel evolution leading tothe necessity of controlled diversion. The sed 
ment transport characteristics of channel reaches in the study area provide : 


2 Old River Diver sion Contr ol a Symposium iT ransactions ASCE V ol 123 19! 
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xcellent example of sediment transport capacity imbalance in an interconnected 
stem of rivers. For the purposes of this paper, the discussion will be con- 
ined to Old River and the upper 41 miles of the Atchafalaya River. 


ATCHAFALAYA RIVER 


The Atchafalaya River flows, for the most part, through the leveed Atcha- 
alaya Basin Floodway which averages about 17 miles in width. At project de- 
ign flood stage this floodway will be required to pass 1,500,000 cfs which is 
0% of the project design flood flow at this latitude. The river channel is in 
he approximate center of the floodway and is leveed from Mile 5 to 53, meas- 
ring from its head. This leveed reach varies from 2,500 ft to 3,500 ft in width 
nd is flanked on the west side by the West Atchafalaya Floodway and on the 
ast Side by the Morganza Floodway. From Mile 53 to 115, flow is channeled 
hrough low-lying swamp and shallow lake areas. 

By virtue of a lesser distance to sea level, complemented by widespread 
low in the lower basin, the water surface slope in the upper leveed reach at 
ppreciable flow is steeper thanthat of the Mississippi River below Old River. 
\lso, the leveed reach is rather deep in comparison with its width. For in- 
tance, in 1957 the bankfull stage width-to-depth ratio of the Atchafalaya River 
t Simmesport was 32 as compared to 120 for the Mississippi River just below 
ne head of Old River. In Old River at Torras, 1.1 miles from the Mississippi 
‘iver, itwas 39. Mean depths measured at bankfull stage in 1957 were approxi- 
aateiy 60, 50, and 45 ft, respectively, at Simmesport, Torras, and Red River 
sanding. 

It cannot be stated as being of exact or uniform occurrence, but there is 
vidence that degradation occurs inthe Atchafalaya River at anytime of medi- 
mand higher stages coincident with low velocities in Old River. At these 
tages Red River sediment supply in the sand range is negligible because of 
lat slopes in the backwater area. Under any circumstance the sand fraction 
f Red River sediment is much finer than that coming from the Mississippi 
rough Old River. Therefore, the important sediment supply which affects 
1e regimen of the Atchafalaya River comes from Old River. Conversely, when 
ere is appreciable flow in the Atchafalaya coming almost entirely from the 
lississippi River, some aggradation of the Atchafalaya is to be expected. 

Extreme cases of the latter condition are infrequent, the last occurrence 
eing in June, 1943. At that time flow in the Atchafalaya was in excess of 
00,000 cfs and all of it was coming from the Mississippi via Old River. Av- 
rage velocities in Old River at that time reachedas high as 8 fps as compared 
ith less than 6 fps in the Atchafalaya. Roughly, this would make the supply 
f sediment to the Atchafalaya about 50% in excess of its transport capacity. 
bservations indicate that at bankful stage the mean bottom elevation at the 
immesport discharge range may raise or lower from 2 ft to 4 ft, depending 
n the percentage of Atchafalaya River discharge being contributed through 
ld River. 

In 1952 a study was made of bed material transport capacity in the Atcha- 
aya River from its head to Krotz Springs (Mile 41). This study was based 
11949-50 channel conditions and water stages observed during the 1950 flood. 
hese data showed that water surface slopes in this reach were not uniform 
ither with respect to stage or location, nor was the channel section uniform 
| average width and depth. In the downstream direction the channel decreased 


oe 


40 June, 1960 HY 6 


in width, increased in depth and, in general, there was an increase in water 
surface Slope. It was indicated that the slope also increased with increasing 
stage. 

76 study the effects which these varying channel characteristics might have 
on sand transport, the river was divided into three reaches as shown in Fig. 2. 
This division of reaches, Barbre Landing to Simmesport, Simmesport to Mel- 
ville, and Melville to Krotz Springs, was dictated by location of gaging stations 
at which daily observations are made. 

Sand transport rating curves were then developed for each of these reaches 
in general accordance with the bed-load function theory developed by H. A. 
Einstein.? At that time the bed-load function theory was relatively new and 
initial studies were directed towards determining its applicability to the Atcha- 
falaya River. Sediment data from frequent observations at the Simmesport 
discharge range were available for this purpose. It was found that application 
of the principles of the bed-load function theory, with a few modifications, would 
produce results in reasonable agreement with the measured sediment loads. 

Fig. 3 shows the sand transport rating curves developed by this method for 
the three reaches. It is interesting to note that the rating curves show the 
leveed reach of channel to be, in general, in a state of equilibrium at flow up 
to 400,000 cfs. This is approximately bankfull flow in the leveed reach. Above 
400,000 cfs, the transport capacity increases in the downstream direction. The 
difference between Reaches 1 and 2 is more or less constant, but in Reach 3 
itis increasingly greater as the discharge increases. Thus, it is indicated that 
at that time there was sediment transport capacity imbalance in the upper 
reaches of the Atchafalaya River at stages above bankfull. 


OLD RIVER 


Old River is approximately 6.5 miles in length and is the connecting channel 
between the Atchafalaya and Mississippi Rivers. Its mean depth at bankful 
stage is less than that of the Atchafalaya andsomewhat greater than that of the 
Mississippi. Formerly its width was several hundred feet less than the Atcha- 
falaya, but in recent years its rate of enlargement has exceeded that of the 
Atchafalaya and now there is little difference between the two, in this respect 
The amount of sand supplied to the Atchafalaya River from outside sources 
includes that which passes into Old River from the Mississippi, plus and minus 
fillor scour in Old River, plus that contributed by the Red and Ouachita Rivers 
system. The latter system is not considered an important contributor of ma- 
terial that will affect the regimen of the Atchafalaya River. 

Discharge through Old River is not entirelya function of stage; itis, rather 
dependent on relative stages at its head and mouth. A fast rising Mississipp 
coincident with small contribution from the Redand Ouachita system produce: 
high velocities inOld River. These velocities are sometimes much higher thai 
those attained in the Atchafalaya at the same time. Conversely, a rapidly fall. 
ing Mississippi following a flood crest results in low velocities in Old Rive: 
with sustained high velocities in the Atchafalaya. Old River flowis also affect 
ed by appreciable flow from the Red and Ouachita system. In 1945a large con 
tribution from the latter system coincident with a falling Mississippi resulte 
in a reverse flow in Old River. This was not uncommon prior to 1945. In ; 


3 U.S. Dept. of Agric., SCS, Tech. Bulletin No. 1026, September, 1950. 
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few instances since that date flow has been reduced to very small amounts, but 
there has been no recurrence of reverse flow. This is attributable to the in- 
creased flow capacity of the Atchafalaya River. 


SEDIMENT TRANSPORT 


As noted previously, the interrelated river reaches under discussion pro- 
vide an excellent example of sediment transport capacity imbalance. To present 
the magnitude of imbalance it is necessary to state sediment loads quanti- 
tatively. Also, in order to show the degree of imbalance for a full range of 
water-type years, a number of years must be considered. This involves la- 
borious computations if carried out in great detail. 

It was not considered necessary for purposes of illustrating this imbalance 
to include the entire bed material load, nor was it considered necessary to 
make detailed load computations. Mechanical analyses of a number of bed- 
material samples taken in the study reach showed medium sand (0.25 to 0.50 
mm) to be the dominant size, comprising approximately 35% of the total. Since 
any grain size group would serve to show relative transport capacities quanti- 
tatively, determinations were made only for the medium sand fraction. For 
analytical load determinations, it was found that a close approximation of the 
results of detailed computations could be attained by simple formulas requiring 
only knowledge of the average velocity and mean depth. These formulas, applied 
to velocities and depths, as measured at frequent intervals at the Torras and 
Simmesport discharge ranges, are the bases for the quantitative medium sand- 
loads presented. 

An examination of mean annual discharges of the principal contributary 
rivers led to the adoption of three water-type years for study. These were 
classified according tothe magnitude of average annual flowin the Mississippi 
River as high-, medium-, and low-water years. For each of these types three 
years were selected as representative of various combinations of mean annual 
flows in the contributary streams. The years 1950, 1945, and 1951 were se- 
lectedas representative of high-water years; 1957, 1958, and 1946as medium- 
water years; and 1944, 1953, and 1954 as low-water years. 

In the high-water years the average daily Mississippi discharge was from 
826,000 cfs (1951) to 879,000 cfs (1950); for medium-water years from 605,000 
cfs (1946); to 685,000 cfs (1957); and for low-water years from 333,000 cfs 
(1954) to 567,000 cfs (1944). Variation in average daily contributary flow from 
the Red and Ouachita Rivers system, that has a large influence on flow in Old 
River, was from 46,000 cfs (1951) to 102,000 cfs (1945) for the high-water 
years; 70,000 cfs (1958) to 88,000 cfs (1957) for the medium-water years; and 
from 25,000 cfs (1954) to 55,000 cfs (1944 and 1953) for the low-water years. 

Although 1944 was placed in the low-water year category, it is a borderline 
case from the standpoint of Mississippi River flow and could just as well have 
been placed in the medium-year category. However, from the standpoint of 
Red-Ouachita and Atchafalaya Rivers flow, it is comparable with 1953, which 
is a true low-water year. A tabulation of average daily discharges and com- 
puted sediment transport capacities for the years selected for study is shown 
in Table 1. : 

Since sediment transport capacity is directly related to stream velocity, a 
comparison of daily average velocities at Simmesport and Torras for the se- 


a . : 
lected years is shown in Figs. 4, 5, and 6. Fig. 4 shows the high-water years 
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1950, 1945, and 1951. In each of these years the annual medium sand transport 
Capacity of the Atchafalaya River exceeded that of Old River. The greatest 
difference was in 1945 and amounted to nearly 4 million tons. The next highest 
was 1950 with nearly 3 million tons, followed by 1951 with one million tons. A 
comparison of the velocity graphs shows why the difference was so much 
greater in 1945 than 1951. 

The volume of water passing Natchez in 1951 was about equal to that in 
1945, but the peak discharge was much higher in1945. This is reflected in the 


TABLE 1,.—AVERAGE DAILY DISCHARGES AND COMPUTED SEDIMENT 
TRANSPORT CAPACITIES 


Medium Sand? 
Transport Ca- 
pacity in 
Million Tons 


Excess Ca- 
pacity of Atch. 
River over 
Old River 
(Medium sand 
fraction in 

Tons) 


Average Daily Discharge in 
efs x 1073 


Missis- Atchafa- 
sippi Atchafa-| Old | Ouachita Old laya 
River | laya R. | River | and Red©| River River 


High Water Years 


86 5.699 
832 264 147 102 3.396 7.233 
46 3.925 


2,832,000 
3,837,000 
1,030,000 


1957 3.527 (3) (4)-(8) 531,000 


4,058(4) 


22 Apr-3 Aug 1.234 2,824,000 
25 Oct-31 Dec 0.463 | 1.168 705,000 
1958 4,094 | 6.192 2,098,000 
1946 2.304 | 4.542 2,238,000 


1944 1,893,000 


1953 - 
i May-31 Jul 1,565,000 
| Aug-31 Dec -61,000 
1954 242,000 


2 Bed material samples show the average medium sand fraction to be about 35%. 
[This percentage was used in determination of transport capacities. 

b Above Old River. 

c principal tributaries to the backwater area. 


hydrographs shown in Fig. 4 for the Atchafalaya River, the 1945 peak being 
630,000 cfs and in 1951, 425,000 cfs. Also in 1945 there were sustained high 
flows from the Red-Ouachita system which was a deterrent to flow through 
Old River. Consequently, Atchafalaya River velocities were greater for 90% 
of the period from January 1 to August 17th. Old River velocities were domi- 
nant for the most partonly during the low-water season at which time sediment 
loads are relatively small. 

_ The year 1951 was characterized by sustained periods of near-bankfull flow 
in the Mississippi River. In 8 of the 12 months its flow equalled or exceeded 
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1,000,000 cfs. Flows less than 300,000 cfs occurred in only 4 days during the 
year. By contrast, in the year 1945, there were only 5 months in which the flow 
equalled or exceeded 1,000,000 cfs, and there were 53 days when the flow was 
less than 300,000 cfs. In the year 1951 there were no sustained periods of ap- 
preciable flow from the Red and Ouachita system. Therefore, flow in Old River 
was generally steady, and its mean stream velocity compared favorable with 
that at Simmesport throughout most of the year. As a result, the excess sedi- 
ment transport capacity of the Atchafalaya River was much less than in 1945. 

The year 1950 was characterized by more than usual dominance of Old River 
velocity during both the spring rise and the low-water season. However, this 
was more than offset by greater velocities in the Atchafalaya River during most 
of the high-water period. The superior velocities and large volume of flow in 
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the Atchafalaya River at this time accounted for over 50% of its 2,832,000 tons 
total excess for the year. 

Fig. 5 shows comparative velocities in the medium-water years 1957, 1958, 
and 1946. In each of these years the annual medium sand transport capacity 
of the Atchafalaya River exceeded that of Old River. In 1957 the higher ve- 
locities in the Atchafalaya River in May, June, and July were sufficient to per- 
mit the Atchafalaya River a sediment-transport-capacity during this period 
alone to exceed that for the entire year in Old River. This illustrates the po- 
tential of circumstances for great imbalance in sediment transport capacity. 
This year is also characterized by the abnormal situation of higher velocities 
in the Atchafalaya River during most of the low-water period. In both 1958 and 
1946 the excess capacity of the Atchafalaya River was more than 2 million tons 
of the medium sand fraction. 
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Fig. 6 shows comparative velocities in the low-water years. As in other 
water-type years, the Atchafalaya River’s capacity for transporting sand ex- 
ceeds that of Old River in each of these cases. Here again the potential for 
great imbalance is illustrated. For example, in May, June, and July, 1953, 
the transport capacity of the Atchafalaya River was nearly 7 times that of Old 
River. Further, in those 3 months the transport capacity of the Atchafalaya 
was over 7 times that of Old River for the 5 months, August through December, 
a low-water period during which velocities in Old River consistently exceeded 
those in the Atchafalaya River. 

In 1944 the excess transport capacity of the Atchafalaya was comparable to 
some medium-water years, and nearly twice that of high-water year 1951. 
Year 1954 was an exceptionally low water year. Even so, and in spite of su- 
perior velocities in Old River during most of the year, the Atchafalaya River 
showed a slightly greater annual capacity for transporting medium sand. The 
balance infavor of the Atchafalaya River for this year was effected bythe May 
Mississippi River rise concurrent with increased flow from the Red and 
Ouachita system. From the first of May to the fourth of June the sediment 
transport capacity of the Atchafalaya River exceeded that of Old River by ap- 
proximately one million tons. From June 15 to October 10 the Old River ca- 
pacity exceeded that of the Atchafalaya by approximately 300,000 tons. This 
shows, again how effectively a short period of higher velocity and discharge in 
the Atchafalaya River can offset long periods of lower discharge which are 
favorable to superiority of transport in Old River. 


CONCLUSIONS 


It is conclusive from the data presented that the distributary system of the 
Mississippi River in the vicinity of the latitude of Old River is in a state of im- 
balance in sediment transport capacity. This imbalance has been found to exist 
within the leveed portion of the Atchafalaya River itself during large floods, 
and more importantly in its relation with supplying streams in all types of 
water years from low to high, when considered on an annual basis. 

The end result of this state of imbalance can lead in only one direction— 
continuous annual enlargement of the Atchafalaya River regardless of water- 
type year. In substantiation of these findings from the standpoint of sediment 
transport capacities, it is a matter of record that the rating curve for the 
Simmesport discharge range, during the period 1951-58, lowered 2 ft at ap- 
proximate bankfull discharge. In this seven-year period there were no high- 


water years; three of the years were of medium-year magnitude and four were 
low-water years. 
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DRAG AND LIFT ON SPHERES WITHIN CYLINDRICAL TUBES 


By Donald F. Young! 


SYNOPSIS 


A method for evaluating the lift and drag forces that are exerted by a mov- 
ng fluid on a spherical particle resting on the bottom of a cylindrical tube is 
escribed. Experimental data are presented to show the variation in the lift 
nd drag force with the particle to pipe-diameter ratio and the Reynolds num- 
er. 


INTRODUCTION 


One of the most important aspects of the general problem of the flow of 
uSpensions is the study of the conditions necessary for maintenance of the 
Jlid particles in suspension. This specific problem is, in turn, related to 
le evaluation of the fluid forces that act on the individual particles. The re- 
iltant dynamic fluid force acting on a particle can be resolved into two com— 
nents, a drag force parallel to the direction of flow, and a lift force normal 
) the direction of flow. 

Although some attention has been given to the effect of the drag component 
1 the behavior of suspensions, little work has beendone in connection with the 
termination of the lift component. It is probably true that in certain prob- 
ms the lift force has little significance, however, for the study of incipient 
otion of particles resting on the bottom of channels or pipes the lift force 


Note.—Discussion open until November 1, 1960. To extend the closing date one month, 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ty of Civil Engineers, Vol. 86, No. HY 6, June, 1960. 

1 Assoc. Prof., Dept. of Theoretical and Applied Mechanics, Iowa State Univ. of 
ience and Tech., Ames, Iowa. 
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may be of primary importance. S. Leliavsky has summarized2 the work that 
has beendone inthis area and points out that it is “an almost unexplored aspect 
of the problem.” The purpose of the investigation reported on herein was to 
develop a method for evaluating the lift and drag forces that act on a single 
spherical particle located on the bottom of a tube, and to study the parameters 
that affect the magnitude of these forces. 


MEASUREMENT OF LIFT AND DRAG FORCE 


The free-body diagram for a spherical particle in contact with a fixed 
boundary inclined at an anglef is shown in Fig.1. The resultant dynamic fluid 
force has been resolved into a lift component L and a drag component P. I 
should be noted that, in general the resultant fluid force does not pass througt 
the center of gravity of the particle. C.M. White has investigated3 the loca- 
tion of the drag component in some detail. The weight of the particle minus 


FIG, 1.—FREE-BODY DIAGRAM 


the buoyant force isdenoted by W'(W’ > 0) and theforces Fand N are the com 
ponents of the boundary reaction on the particle. 
Equilibrium equations for the particle are 


L= W' cos. p=N >... « s aki see (1 

in the n - direction, 
P=W' sinp- F i ike cents ob Sa 
in the t - direction, and . 
° s 
PLAT Fic c eee cee eee eee eee al 


wee Introduction to Fluivial Hydraulics,” by S. Leliavsky, Constable and Co., 195 
p. 64. 
3 “The Equilibrium of Grains on the Bed of a St x I i i 7 
ream,” by C. M. White, Proceed 
Royal Soc. of London, A174, 1940, p. 322. : 7 
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for moments with respect to an axis through the center of gravity of the par- 
ticle of diameter d. 

By varying the rate of flow past the particle, the drag component P can be 
changed and the particle made to move up or down the plane with a constant 
velocity or to remain stationary. As the angle Bis increased, a greater ve— 
locity is required to hold the particle stationary or move it up the plane. How- 
ever, as 8 increases, the component of the apparent weight normal to the pipe 
wall W',, that opposes the lift force decreases. Therefore at some critical 
angle it should be possible to have the lift force equal W"y so that any slight 
increase in the value of L will cause the particle to rise off the bottom of the 
pipe. Since W‘y can be readily measured, the lift force can be determined for 
this condition of impending motion. Alsofor this condition of impending motion 
the drag force is equal to W"sin Bsince F = 0. Thus both the lift and-the drag 
force can be determined. This method of determining these forces has the ad- 
vantage that no attachments need be made to the particle to hold it in position. 

However, it should be noted that as motion impends in the direction of the 
lift force both the normal force N and the frictional force F approach zero. 
Therefore in order that no rotation take place, Eq. 3 indicates that the product 
P Z must approach zero. Since P + 0 this implies that either Z = 0 or that Z 
approaches zero as the lift force approaches W’,. It is thus questionable 
whether impending motion in the direction of the lift force will occur without 
rotation. However, actual observations indicated that a particle can be made 
to rise off the bottom without rotation indicating that the necessary conditions 
were satisfied, at least for the tests run in this investigation. 


EXPERIMENTAL PROCEDURES 


The experimental set-up used to measure the lift and drag forces by the 
procedure outlined was quite simple. Itconsisted primarily of a constant-head 
tank connected to a glass tube which could be inclined at various angles. The 
angles could be measured to within 1.0°. Parallel needle valves were used to 
2ontrol the rate of-flow through the system. The fluid temperature remained 
rélatively constant for all tests, and the fluid used in all tests was water. A 
schematic diagram of the equipment is shown in Fig. 2. 

It was expected that for turbulent flow in the pipe the motion of the particle 
would be erratic and measurements difficult to obtain. For this reason it was 
lesired to run all tests under laminar-flow conditions in the pipe. Therefore 
1 particle of low specific weight was needed so that the desired conditions 
“Ould be obtained with relatively low velocities. To meet this requirement 
ylon spheres with a specific weight of 1.1923 g per cucm were used. The 
est facility was arranged so that various sizes of glass tubing could be used. 
n order that the flow at the measurement station be fully developed laminar 
low, this station was located at a distance of approximately 4 ft from the en- 
rance and 2 ft from the discharge end of the tube. A commonly used equation 
or estimating the transition length x is 


SF O2B tas Raney ats 7 ewes er a4) 
*, 
n which r is the tube radius and R is the Reynolds number based on the tube 
adius.4 The maximum transition length that occurred during the tests, as 


ee a ee REAR Colby Wi. IGM a a cies. otk bens, 
_ 4 «Applied Hydro- and Aeromechanics,” by L. Prandtl and O.G. Tietjens, Dover Pub- 
ications, Inc., New York, 1957, p. 22. 
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estimated from Eq. 4, was 3.3 ft. It was therefore assumed that for all the 
tests the flow was fully developed at the measurement station. Table 1 sum- 
marizes the combination of particle and tube sizes used in the experimental 
investigation. 

For a given test the angle B was fixed and the motion of the particle in the 
tube was controlled by the rate of flow from the constant-head tank. As ex- 
pected, atsome critical angle the particle could be made to rise off the bottom 
of the tube with essentially no motion parallel to the longitudinal axis of the 
tube. The particle would lift off the bottom, then the lift force would apparent- 
ly decrease, and the particle would drop back so that there was actually a “hop- 
ping” motion. 

The procedure followed in each test was to adjust the angle and rate of flow 
until the particle had impending motion normal to the pipe wall and to record 
the data for this condition. The rate of flow in the tube was obtained by col- 
lecting the fluid as it discharged from the tube and determining the time re- 
quired to collect a given volume. 


FLEXIBLE 
TUBING 
To 
DISCHARGE 
TANK 
FLEXIBLE 
TUBING 


FIG, 2.—EXPERIMENTAL APPARATUS 


One disadvantage to this method of measuring the fluid forces acting on 
particle is that a certain amount of personal judgment is involved in establish 
ing the condition of impending motion of the particle. In order to have som 
index to the reliability of the resulting data, each test with a particular com 
bination of particle and pipe size was repeated at least seven times and th 
standard deviation of the lift and drag coefficients was evaluated. The stand 
ard deviation from the mean values of the lift and drag coefficients was les 
than 3% of the mean values for all tests. As indicated in Table 1 the range « 
values for the (d/D)-ratio was from 0.282 to 0.635. The Reynolds numbe 
range, based on the pipe diameter and mean velocity, was 360 to 1,117. 
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TABLE 1.—SUMMARY OF EXPERIMENTAL RESULTS 


—— + a = ie r 

da d, in Dan 

D em em Rq Rp Cp Cy, 

] + 

282 0.321 1.135 220 780 1.50 0.61 
345 0.480 1.389 385 1,116 1.61 0.73 
O53 0.401 1.135 284 805 WES 0.78 
oa 99 0.321 0.803 166 416 2.58 “leas Kes 
403 0.559 1.389 450 i Hes Be Ie 1.86 0.87 
1.423 0.480 L135 330 780 Dit t: 1.04 
1.459 0.638 1.389 495 1,078 2.26 1.10 
493 0.559 1.135 380 TAL 2,58 wet 
500 0.401 0.803 204 408 3.30 1.62 
Zo19 OLT2t 1.389 538 1,037 2.76 1.36 
562 0.638 a ESS 402 715 3.39 Ce 
098 0.480 0.803 215 360 5,02 2.58 
.635 0.721 1.135 411 647 4.67 2.43 


ANALYSIS OF EXPERIMENTAL DATA 


It was assumed that the variables influencing the lift force are the mean 
elocity V, fluid density p, fluid viscosity ., particle diameter d, and pipe di- 
meter D, so that 


Pree (View hd, Des hi ds alg aos (5) 


'q. 5 can be written in terms of three dimensionless parameters as 


here A is the cross-sectional area of the particle and the constant 2 has been 
rbitrarily introduced to put the equation in a more conventional form. There- 
ore the expression for the lift force is 


2 
+, Vv 
L=C, pAyz.-::-: mad tuk sate: OE OCHO erie ta Es aus Ml) 


here the lift coefficient Cy, is given by 


he Reynolds number Rp is based on the pipe diameter. Of course a Reynolds 
umber based on the particle diameter could have been used. However it was 
It that the velocity distribution preceding the particle, which is a function of 
1e pipe Reynolds number, would have considerable influence on the lift and 
rag coefficients and thus Rp was used. However the results can readily be 
iterpreted in terms of the Reynolds number Rg based on particle diameter 


Ince 
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Similarly the drag force can be expressed as 


2 


p Vv 
P =r p Maa ita da ore nse eae (10) 


where 


Ch = %o (r,, A Baad oaschis xn 5m ROeaias wae uiee ante (11) 


Since the lift and drag coefficients can be evaluated from the experimental 
data, the problem reduces to one of determining the unknown functions $1 an¢ 
$9. Table 1 contains all the data and Fig. 3 is a plot of the experimentally de- 
termined values of CI, and Cp versus Rp for the various diameter ratios. Due 
to the limited number of particle and tube sizes available and the manner ir 
which the data were obtained it was not possible to hold either d/D or Rp con- 
stant for a series of tests. It was therefore difficult to correlate the data 
However an examination of Fig. 3 suggested that both the lift and drag coeffi- 
cients were inversely proportional to Rp! 3, It was further assumed that the 
equations for the lift and drag coefficients could be written as 


Ewha d 
Cy => R51/s f, (3) € Poti S -* eae: wee eo eke (12) 
and 
aaelay: d 
Ch = Rpl/3 fy (5) ete C hp) 96 CARE ee ie ee (13) 


It is apparent that if these assumed forms are correct then a plotof Cr Rp! / 
and Cp Rp!/ 3 versus d/D should correlate all the data, and these plots were 
used to indicate the validity of Eqs. 12 and 13. It was found that for the range 
of variables covered in the investigation the data do appear to correlate satis- 
factorily as shown in Figs. 4 and 5. The empirical equations for the drag anc 
lift coefficients which fit the data shown in these figures are 


and 


2.28 aa 
c= xy [9 5(2) -2.50| Pll gOS areal 


It should be emphasized that although the relatively simple form for corre: 
lating thedata appears to be reasonably satisfactory, it is not likely that it wil 
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> valid over the entire range of Reynolds numbers that is commonly of inter- 
st. In this investigation the range of Reynolds numbers based on particle di- 
neter Rq was 166 to 538. It is well knownfrom conventional drag-coefficient 
sts for spheres in an essentially unbounded fluid that this range of Reynolds 
imbers represents a transition region and that a simple valid expression for 
ie drag coefficient in this range cannot be extended into the regions where 
scous forces are predominant or into regions where viscous forces are neg- 
gible. 
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: FIG. 3.—VARIATION OF COEFFICIENTS WITH REYNOLDS NUMBER 


: Figs. 4 and 5 indicate that for a given Reynolds number both the lift and 
ag forces increase with increasing values of d/D and for a given value of d/D 
th the lift and drag forces decrease as the Reynolds number increases. 
ese results for the drag coefficients are in qualitative agreement with those 
tained® by J. S. McNown and J. T. Newlin in their study of the drag force 
ting on spheres located in the center of a cylindrical tube. Data obtained 


POE ES at ae eR 
5 «Drag of Spheres within Cylindrical Boundaries,” by J. S. McNown and J. T. New- 
Bei cccsings, First U. S. Natl. Congress of Applied Mechanics, 1951, p. 801. 
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FIG. 5.—CORRELATION OF LIFT DATA 
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from curves given by McNown and Newlin are plotted in Figs.6 and 7. Itis 
noted that the drag coefficients for spheres located near the wall of a tube in 
which the velocity of approach has a parabolic distribution are considerably 
higher than those for spheres located in the center of the tube in which the 
velocity of approach is uniform. J. Happel and B. J. Byrne derived6 analytical 
expressions for the drag force acting on a sphere in the Stokes range placed in 
the center of a cylindrical tube for both the case with a uniform approach ve- 
locity and with an approach velocity with a parabolic distribution. They found 
that the drag coefficients for the case with the parabolic distribution were 
roughly twice those for the case with the uniform distribution if the character- 
istic velocity is taken as the mean velocity. This difference in drag coefficients 
is of the same order of magnitude as those indicated in Figs. 6 and 7. 


3.0 


CALCULATED FROM 
EXPERIMENTAL DATA 


2.0 
2D 
1.0 
cNOWN & NEWLIN 
0.8 
0.6 
107 2 4 6 8 10 
Rq 


FIG, 6.—VARIATION OF Cp WITH R, FOR D/d = 0.4 


From Figs. 4 and 5 it is apparent that the ratio of the drag to the lift force 
is a function of only d/D. The plot of this ratio versus d/D is shown in Fig. 8. 


CONCLUSIONS 


In this study a method was developed to determine the lift and drag forces 
acting on a spherical particle resting on the bottom of a circular tube. A num- 
ber of tests was performed to confirm the feasibility of using this method and 
to determine the values of the lift and drag forces. It is concluded that the 


CE aa 
6 “Motion of a Sphere and Fluid in a Cylindrical Tube,” by J. Happel and B. J. Bryne, 
ndustrial and Engineering Chemistry, Vol. 46, 1954, p. Inet 
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FIG, 7.—VARIATION OF Cp WITH D/d FOR Rg = 300 


FIG, 8.—VARIATION OF DRAG TO LIFT RATIO WITH d/D 
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proposed technique can be used to measure these forces accurately. In addi- 
tion, Fig. 8 shows that the lift force is of the same order of magnitude as the 
drag force and that the ratio of drag to lift, although a function of thed/D-ratio, 
iS approximately equal to 2. 

The range of values of the parameters studied in this investigation is limited, 
but the same type of equipment and technique can be used to extend the range 
with the proper choice of fluids and particles. Probably the most important 
result of this study is the establishment of the fact that for spheres near the 
wall of a tube alift force with a magnitude comparable to that of the drag force 
does actually exist. It is thus apparent that the lift force should not be over- 
looked in studies related to the incipient motion of particles resting on stream 
deds or pipe walls. 
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MODELS PRIMARILY DEPENDENT ON THE REYNOLDS NUMBER 


By W. P. Simmons, Jr.,1 M. ASCE 


SYNOPSIS 


The design, construction, and operation of models of closed-conduit fluid 
systems, and interpretation of test data are discussed. Viscous and inertia 
forces in the flowing fluids of such models, and hence Reynolds number, are 
oredominant factors. General rules for model size, construction, instrumen- 
ation, operating Reynolds number ranges, and data analysis are given. 


INTRODUCTION 


Hydraulic studies of models of closed-conduit flow systems provide a re- 
iable means for obtaining useful and accurate design data. However, there 
ure pitfalls in the use of models that must be recognized and avoided if the 
lecessary accuracy and reliability are to be obtained. Thus, a knowledge of 
he similitude relatienships between models and their full-sized counterparts 
s essential. An excellent discussion of these similitude factors has been pre- 
sented by H. K. Liu and M. L. Albertson,2 and by H. Rouse;3 therefore, a full 
.ecount of the details will not be given here. It will be in order, however, to 
lefine what is meant by closed-conduit flow and to state briefly the fundamental 
nodel relationships that must be observed. 

Note.—Discussion open until November 1, 1960. To extend the closing date one month, 
. written request must be filed with the Executive Secretary, ASCE, This paper is part | 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 86, No. HY 6, June, 1960. 

1 Hydr, Engr., Div. of Engrg. Labs., Bur. of Reclamation, Denver, Colo. 

2 “Significance and Application of Froude and Reynolds Numbers as Criteria for 
‘imilitude,” by H. K. Liu and M. L. Albertson, presneted at the June 1959 ASCE Con- 


ention in Fort Collins, Colo. 
3 “Engineering Hydraulics,” by H. Rouse, John Wiley and Sons, Inc., New York, 1950, 
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Closed-conduit flow is definedas flow in an enclosed system where no free 
water surface exists, or flow around a deeply submerged body where the effect 
of a free water surface is negligible. Typical examples of structures invoiving 
closed-conduit flow include piping systems, pressure tunnels, control gates, 
pumps, and turbines. 

In studies of these systems, gravitational forces are not a factor affecting 
the flow because the flow direction always will be down the slope of the energy 
grade line and independent of the slope of the pipe itself. Surface wave motion 
and surface tension are nonexistent because there is no free surface, and, 
therefore, are not a consideration. Elastic compression of the fluid, which in 
the case of gases can be significant, is usually negligible with liquid flow. Thus, 
ofall the major forces and fluid characteristics affecting fluid flow, viscosity 
remains as the predominant factor in the usual closed conduit, hydraulic flow 
problem. 
ae or vp is the ratio of the fluid inertia forces 
to the fluid viscous forces. The number is dimensionless, and p is the fluid 
density, V denotes the flow velocity, D is a characteristic dimension (usually 
the diameter in pipes), and y and v are the absolute and kinematic fluid vis- 
cosities, respectively. Values for v for fluids at normally encountered tem- 
peratures are given in Fig. 1. 

Usually the Reynolds number of a prototype structure is very large. Dupli- 
cating this number on scale models is impractical, particularly if the same 
fluid is usedon both structures, for, as the size decreases, velocities must in- 
crease proportionately. With high initial prototype velocities, the problem 
rapidly would seem to become impossible. 

Study of the Reynolds number equation shows that the effects of viscosity 
are particularly important when flow velocities and, hence, inertia forces, are 
low. But when flow velocities are high, inertial forces assume such great sig- 
nificance that viscosity becomes relatively minor. This is a basic factor that 
makes it possible to achieve reasonable similarity between a model and its 
larger-sized prototype. By making the model large enough, or by using ve- 
locities high enough, or by both, inertial forces can be made to predominate 
over viscous forces and a reasonable similarity with the prototype will exist. 
The other requirement is that the model be large enough to insure flow in the 
fully turbulent range whenever full turbulence occurs in the prototype. Thus, 
the procedure followed in making usual model studies of elements of closed- 
conduit hydraulic systems is to build reasonably large, very accurate models, 
and to operate them with flow velocities equal to or higher than scaled veloci- 
ties. Full turbulence is generally encountered on prototype structures, and the 
discussion in this paper will be limited to these conditions. 


The Reynolds number, 


MODEL DESIGN 


A first consideration in designing a model is the selection of the testing 
fluid. Water, the most generally available and perhaps most satisfactory li- 
quid, is commonly used. Water also is usually the fluid used in prototype 
structures and naturally associated with model studies. In many types of closed- 
conduit testing, air can be used as the testing fluid with a real savings in time 
and costs.4 For this reason, low-velocity air tests should be considered in the 


ee ee Ee ee 
4 “Model Tests Using Low Velocity Air,” by J. W. Ball, Transactions, ASCE, 1952 
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planning stages of any model study. Sometimes it may be desirable to take 
advantage of both testing fluids by making rapid, low-cost air tests in the early 
stages of the study, and using water for the final tests and calibrations. 

The model should be made as large as practicable, considering construc- 
tion and operating costs and the benefits that may be derived. Increasing the 
size of a model usually enhances its usefulness and improves its accuracy. At 
some point, however, the increase in cost and the difficulty of operation will 
offset any size advantage. Current practice is to follow precedent when avail- 
able and to erron the side of largeness insofar as space and water supply per- 
mit. In general, most closed-conduit models will fall within the scale ratio 
range of 1:5 to 1:30. Undistorted scale modeling is the rule. 

For reasonable similitude between model and prototype, it is usually neces - 
sary to operate the model at Reynolds numbers of approximately 1 x 106 or 
more. At these values the frictional coefficient f in the Darcy-Weisback for - 
mula 


is constant for many of the most frequently encountered pipe roughnesses, and 
nearly constant for very smooth pipes (Fig. 2). In the case of valve or conduit 
models using scaled heads and velocities, this requires about a 6-in. pipe di- 
ameter. Consideration must be given to the pump capacities to determine whe- 
ther adequate discharge and head will be available for the scale ratios being 
considered. In this regard, it must be remembered that part of the tests may 
be made at reduced gate openings and that heads much larger than scaled may 
be necessary to maintain reasonable values of Reynolds number. 

The extent of the total prototype structure to be represented in the model 
must be determined. Usually, interest is centered around a particular part of 
a system, rather than on the system as a whole. In such cases, it is necessary 
to model only that particular portion of the structure, with enough inlet and 
outlet piping to produce flow conditions similar to those in the prototype. For 
example, in the case of a control gate at the midpoint of a long conduitin a dam 
approximately 20 diameters of model conduit upstream from the gate, the gate 
itself, and about 5 diameters of conduit downstream might suffice (Fig. 3). If 
the flow distribution were severely distorted by multiple bends or a partially 
opened valve at the start of the conduit, a longer upstream conduit would be 
needed. As another example, several control gates might be connected to a 
branching manifold (Fig. 4). In this instance, the characteristics of the flow 
reaching each gate would be affected by the configuration of the manifold and 
the number and location of gates operating. Thus, it would be necessary to 
model the entire manifold and all the gates to represent fully conditions that 
could occur in the prototype. If studies were made of a high-head conduit in- 
let, it might be necessary to provide a pressure tank to represent the reser- 
voir (Fig. 5). In studies with lower heads, a head box may be adequate (Fig. 6). 
Each structure to be modeled must be examined, and engineering judgment 
used, to determine the model limits. : 

Frequently, it will be advantageous to model a critical area only within a 
structure and to do thison an enlarged scale. An example would be a gate slot 
for restraining the leaf of a control gate. A relatively large slot can be formed 
in the wall of a test conduit, and accurate tests can be made if care is taken to 
have the opposite wall far enough away so that there is no hydraulic interfer- 
ence (Fig. 7). Another example would be the use of pie-shaped sectors of cir- 
cular control valves (Fig. 8). Sectors of less than 180° are the most common. 
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Time and unnecessary expense may be saved by selecting the scale rati 
to accommodate standard sized pipes and other model components that are al 
ready on hand. Odd scale ratios, such as 1:17.35 willin no way affect theaccu 
racy of the model or complicate the data analysis. 

Materials ordinarily usedin water models include bronze castings, plate an 
sheet brass, galvanized sheet steel, plastics, concrete, wax, and modeling clay 
In air models, wood and plaster are also satisfactory. When moderately hig 
heads and velocities occur in hydraulic models, rigid construction must b 
used. Models circular in cross section, such as needle valves, use bronze cast 
ings to advantage. When the shapes are rectangular in cross section, machine 


FIG, 9.—TRANSPARENT PLASTIC MODEL 


brass plates, doweled and bolted together, work well. An advantage of th 
build-up design is that model changes can be made readily. Transparent plast 
is used frequently when the flow within the structure must be seen. Mode 
entirely transparent are fairly common (Fig. 9). In any case, the materia 
selected in the design stages should combine durability, dimensional stabilit 
ease of construction, ease of modification, and corrosion resistance. 

The inclusion of measuring equipment such as piezometers is an importa 
design consideration. Often the requirement of getting the instruments in pla 
within the model dictates the type of construction used. Piezometers are rel 
tively cheap when included in the initial model construction, and should be us: 
generously to offset oversights in defining critical areas. Thought shou 
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verefore be given to the necessary instrumentation and how it will be included 
1 the model so the desired information will be obtained. 

The testing fluid must be supplied to the model in a manner representative 
f prototype conditions. Appropriate conduit length must be provided between 
ny Supply system bends or controls, and the model itself, so that theflow dis- 
‘ibution will be developed to the proper extent by the time it reaches the test 
ection. Straightening vanes are useful in reducing the tendency for spiralling 
ow within the pipes. 

The testing fluid may be supplied directly from pumps, or from surge or 
onstant head tanks. Ideally, it would be desirable to have each model supplied 
rough a constant-head tank, but in many cases this is unnecessary and excel- 
mt results are obtained with the flow supplied directly from pumps. The 
1aracteristics of the pumps being used have a considerable bearing on whether 
¢ not a constant-head tank is necessary. If the pumps are used directly, it is 
2st to run them at appreciable heads and to throttle the flow to the model. In 
is manner, better flow stability is obtained. 

The waste from hydraulic models must either be returned to laboratory 
sservoirs or to the rivers or lakes from which it was drawn. This may be 
ne in the most convenient manner and has little effect on the model operation 
proper outlet conditions are maintainedat the model test section. Discharg- 
3 from air models usually present no problems. 

Accurate and complete drawings are essential for expediting model con- 
ruction and preventing time consuming, expensive errors. Drawings should 
mtain sufficient informationto enable craftsmen to build a structure that con- 
rms to the designers’ specifications. Special drawings of details are often 
-cessary. All informationfor the installationof instruments suchas piezome- 
rs, pressure cells, velocity measuring devices, etc., should be included, 
ereby obviating extensive changes after the model is completed. In fortunate 
ses where craftsmenare experienced in hydraulic model construction, many 
ructural details may be omitted and much of the construction method left to 
e shop foreman. 


MODEL CONSTRUCTION 


Extreme constructionaccuracy is essential for obtaining goodresults from 
ly hydraulic or air model. Care must be taken that the materials, methods, 
id craftsmenare satisfactory and that tolerance checks are made frequently. 
reatestaccuracy should be maintained where abrupt changes in flow direction 
‘cur, and where flow velocities are high or pressures are low. The stage of 
e testing program affects the type and accuracy of construction. In early 
sts, during which many schemes may be tried to determine feasibility, con- 
ruction need not be as careful as in the final stages, when data are used to 
scribe prototype operating characteristics. 

In hydraulic models of high head gates and valves, relatively thick metal or 
astic sections are usedand all important flow surfaces are machined. Thus, 
se tolerances and accurate component alinement can be obtained. In air 
»dels, wood, sheet metal, wax, and lighter plastic sections are permissible 
cause the internal forces and the weight that the models sustain are small. 
ie need for accuracy remains unchanged. For most water models made by 
> Bureau of Reclamation, United States Department of the Interior (USBR), 
insitions from one shape or area to another, manifolds, bends, etc., are 
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formed from 16 to 22 gage sheet metal, or 1/10 in.to 1 /4 in. thick, hot-worke: 
sheet plastic. In cases where flow conditions within the system are to be ob 
served, the transparent plastic sections are ideal. They are usually compara 
tively expensive because wooden patterns must often be made first, and th 
heated plastic formed over them, cooled, trimmed, joined, and flanged int 
final condition. The handwork increases costs, and the advantage of visual in 
spection must be balanced against this extra expense. 

The inclusion of piezometers often complicates the model construction an 
may restrict the choice of materials. In the case of castings or heavy meta 
sections, piezometer holes may be drilled through the bodies and connectin 
taps can be threaded or soldered into the exterior surfaces. When a number © 
piezometers are grouped within a limited area, it may be necessary touse de 
vious passages to prevent one piezometer lead from interfering with another 
In sheet metal construction, metal tubing is inserted through the metal walls 
soldered in place with the tubing normal to the wall, and then smoothed flus 
with the flow surface. All burrs must be removed from the interior of th 
piezometer opening.» When concrete, plaster, wood, or wax is used for form 
ing intricate flow surfaces, it is desirable to embedin the material metal tube 
that extend to the surface, so that the piezometers will have durable, smoott 
square-edged openings. When materials are suitable and superior craftsman 
ship is available so that chipping or damaging donot occur, holes can be drille 
directly into the material. Usually the piezometers should be 1/16 in. or large 
in inside diameter and should extend at least 2 diameters from the openin 
without change incross section or alinement. Piezometers smaller than 1/1€ 
in. diameter may be used, and may be desirable on curving surfaces or wher 
many are crowded into a small area, but they are prone to becoming plugge 
with small bits of material, and tend to induce damping in water models th: 
reduces the sensitivity of measurement. The smaller piezometers work we 
in air models. 

Flexible tubing about 1/4-in. inside diameter is usually used to connect tl 
piezometer taps to the pressure measuring instruments, usually water mal 
ometers. Rubber tubing has been used extensively, but is relatively short-live 
and has anopaque quality that prohibits observing interior conditions. Recen 
ly, flexible, transparent plastic tubings have become available and are high 
recommended. An advantage of the transparent tubing is that entrapped a 
bubbles can be detected easily and removed from the line. 

When measurements are to be made of rapidly fluctuating pressures, ca 
must be taken in matching piezometer diameter, pressure lead diameter, pre: 
sure lead length, the pressure cell itself, and the recorder.9,7 Failure to | 
this may result in obtaining records that reflect serious distortions of the a 
tual pressure history. Deductions made in all innocence and confidence fre 
such records would, of course, be erroneous and misleading. In general, 
should be remembered that the pressure leads should be as short as possib! 
rigid, and of at least 1 1/6-in. inside diameter. All entrapped air and g 
bubbles should be bled from the system, and the temperature of the cell shot 


: . Sb Investigation,” by C.M. Allen and L. J. Hooper, Transactions, ASI 
yd 54-1. 
6 “A Note on the Evaluation of Design of Transducer for the Measurement | 
Dynamic Pressures in Liquid Systems,” by J. R. Barton, Instrument Notes, Stath 
Labs., No. 27, June, 1954, . 

7 “attenuation of Oscillatory Pressures in Instrument Lines,” by A. S. Iberall, | 
search Paper RP 2115, Natl. Bur. of Standards, Vol. 45, July, 1950. 
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Provisions should be made for installing any other measuring equipment 

ded on the model. For instance, it may be desirable to make velocity tra- 
rses across the conduits, and pitot tubes or hot wire anemometers may have 
be inserted through the walls. In other cases, the injection of dyes or fine 
reams of air bubbles may be used to show flow patterns. 

Flow-measuring devices such as Venturi meters, orifice meters, or weir 
mes are required for determining the rate of flow. Much depends on the ac- 
racy of these flow measuring devices and carefully calibrated, properly placed 
uipment is essential. Quality metering devices are available and the se- 
ctions can be based on accuracy desired, cost, convenience of use, and ex- 
cted service life. 

Properly fitted, tightly connected pipes, fittings, and transitions should be 
ed between the modeland the flow-measuring device to insure that all mea- 
red flow reaches the model test section. This requirement applies equally 
water and air models. Leaks in air models are often difficult to detect be- 
use they leave no visible trace. In water models, outside leakage and drips 
e easily detected, but leakage of air into zones of subatmospheric pressure 
less obvious. 

In assembling the model components accurate alinement must be maintained. 
Jequate supports, provided beneati the model, should be placed to provide 
inimum interference when alterations to the model are made. Usually it is 
mvenient to place the model at waist height for ease of construction, oper- 
ion, and observation. 


MODEL OPERATION 


The testing program must be carefully planned to make full use of the model. 
‘eliminary runs are necessary to reveal defects, leaks, or other inadequacies. 
rst tests also indicate whether or not minor changes are needed to make the 
odel and the instrumentation better fitted for obtaining the desired data. This 
ase should not be hurried because it is economically important to be certain 
e model really represents prototype conditions and that the instrumentation 
satisfactory. It is during this adjustment phase that the experimenter be- 
mes acquainted with the peculiarities of the model andadept in its operation. 
ssting should consist of examining systematically each proposed, feasible 
sign for possibilities for improvement in flow conditions, cost reduction, and 
se of maintenance. The experimenter must exercise patience, imagination, 
d ingenuity and be capable of interpreting the model results correctly. The 
ta should be analyzed concurrently with the testing to prevent accumulation 
unnecessary data, and to permit intelligent decisions concerning succeeding 
eps in the testing program. A relatively complete and detailed diary of the 
udy is important. This diary should include dates, general notes on perfor- 
ance, changes made, the model data, and photographic records. Often an ap- 
eciable period of time will elapse between the time of a model study and the 
ne when a report is prepared. Without a record such as the diary provides, 
any factors could be forgotten or overlooked. 

Measurements of pressures, discharges, velocities, etc., must be made 
th precision. Much depends on the values recorded on the data sheets, and 
on pertinent remarks such as pressures being particularly unsteady at cer- 
in locations. The type of pressure and velocity records made depends upon 
2 use intended for the data. In many cases, small pressure fluctuations are 


—_ 
ee 


70 June, 1960 HY 


not significant and the average pressures are the ones desired. Damping de 
vices may be used in these pressure lines to make the reading simpler. Othe 
cases willoccur where transient pressures are desired. These measurement 
must be made with special equipment such as pressure cells and electroni 
recorders that are capable of following the pressure swings. No appreciabl 
damping should occur during these measurements. In all pressure readings 
datum points or gage zeroes must be recorded so that actual pressure value 
at the piezometers will be known. The testing must be extended over the en 
tire range of equivalent discharges and operating conditions that will occur j 
the prototype. Frequently it is desirable to extend the range of the tests t 
heads much higher than the equivalent heads anticipated for the prototype t 
fully exploit the model. 

In addition to routine testing of immediate interest, the experimenter shoul 
be alert to obtain general information that will be of value at a later date. Re 
sults of experimentson many models are required for assembling general de 
sign information. Time and money may be wasted if all the useful informatic 
is not obtained from each model under test. 

Studies made using air as a flowing fluid usually require that the flow ve 
locities be maintained below approximately 300 fps. If higher velocities a1 
used, compressibility would have to be considered. The density of the air wi 
vary from day to day with changes in temperature and/or barometric pressur 
Thus, if pressures are read on water manometers in feet, and are converte 
to feet of flowing fluid (air), the density of air relative to that of water must |} 
known. The air temperature and barometric pressure must therefore be re 
corded during air studies. ‘ 


INTERPRETATION OF RESULTS 


If care is taken to operate the model at Reynolds numbers of about 1 x1 
or more, the data may be analyzedusing the same relations as for models bas 
on the Froude law. Thus, direct geometric scaling of linear values is possibl 
If the model were run with scaled heads and discharges, the model pressur 
would be multiplied by the scale ratio to obtain prototype pressures. Tit 
intervals and flow velocities are obtained by multiplying model values by t 
square root of the scale ratio. The discharge is obtained by multiplying t 
model flow by the five-halves power of the scale ratio. 

If the model is tested at heads and discharges that differ from scaled prot 
type values, the data may be placed in a dimensionless form and then appli 
to prototype conditions. An example is shown in Fig. 10. A typical dimensio 
less pressure factor relation is 


(hx = ho)/ Vo"/2 8) 


where hy is the head at the piezometer in question, hg is the head at a refe 
ence piezometer, and Vo is the Q/A velocity at the reference station. 1 
reference station is usually taken at a convenient point just upstream from 
test section itself. Other pressure factor relations may be used as conditic 
dictate, provided they involve a pressure relationship from a base point an 
parameter, such as the velocity head, that describes the flow. 
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Data in the foregoing dimensionless form is applied as follows: Determine 
from the relation 


the rate of flow Q in the prototype structure. In Eq. 2 Ais a significant area 
in the gate or conduit, Cg denotes the coefficient of discharge, and H is the 
total head differential across the structure. After obtaining Q, compute the 
velocity head at the reference station. The pressure factor value for the par- 
ticular piezometer is then multiplied by this velocity head to give the pressure 
change in feet of water from the piezometer to the reference station. The pres- 
sure at the reference station is found by adding any back pressure on the struc- 
ture to the head caused by the presence of the structure, and given by the re- 


v2 
lation (~ +t) = The computed pressure change, added algebraically to 
Cq 8. 

the reference pressure, gives the prototype piezometric pressure to be expect- 
ed. Data obtained from air model tests are best expressed in pressure factor 
form, using care that all values are expressed in feet of the flowing fluid. The 
pressure factors are dimensionless and are applicable directly to prototype 
conditions using other fluids. 

Often negative, or subatmospheric, pressure heads are obtained from models 
which, when scaled directly, yield values too low to be possible in the proto- 
type. An example would be a 3-ft subatmospheric pressure head measured on 
a 1:20 scale model. Direct scaling yields a prototype head of 60 ft below at- 
mospheric pressure. This is unrealistic because no pressures lower than ab- 
solute zero can be obtained. At sea level, this is about -34 ftof water. Actually, 
water begins to boil when the pressure is reduced to vapor pressure, about 0.7 
ft absolute, and cavitation occurs. Thus, any scaled-up model pressures thai 
are subatmospheric to the extent of about 30 ft to 33 ft indicate that cavitatior 
willoccur in the prototype. If cavitation pressures are indicated over appreci- 
able areas of the structure, large-scale separation is likely and the flow pat- 
tern may be materially changed. Specialized testing is necessary to evaluate 
such structures. In the usual case, the mere suggestion of pressures near th 
cavitation range is sufficient to justify redesign so that noise, vibration, ant 
damage caused by cavitation will not occur. 

In some cases, negative pressures of limited extent may be acceptable i 
prototype structures. When flow velocities and pressure fluctuations are moder 
ate or small, indicated pressures 15 ft below atmospheric are not considere 
objectionable. It is usually dangerous to allow pressures lower than this be 
cause unexpected surface roughnesses, vorticity, and flow turbulences ma 
momentarily lower local pressures to the cavitation range. 

Studies of cavitation potential on particular structures are frequently mad 
and often require special test equipment and instrumentation because of th 
complexity of the problems. Certain parameters are commonly used to expres 
cavitation potential on a structure. One is the cavitation index 
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n which ho is the pressure at a reference point, h,, is the vapor pressure 
elative to the atmosphere (negative), and Vo is the velocity at a reference 
oint. Another index especially useful for systems such as valves in pipelines 
=" 


SES ge oe ng A eS Oe eee (4) 
Hy Sh, 


nwhich ho is the pressure in the pipeline well downstream (for instance adis- 
ance of 12 D), hyp is the vapor pressure, and HT is the total head upstream. 
‘his parameter has the advantage of describing pressure conditions which 
aust exist with given flow velocities and valve characteristics if cavitation is 
> be avoided. With either parameter, for the same value of K in model and 
rototype, the pattern of cavitation will be the same. But since frequency of 
ibration and rate of pitting vary with velocity, care must be used in predict- 
ag prototype structural behavior.9 

The head loss through the system is obtained by taking pressure and velocity 
1easurements at appropriate points in the system. Similarly, the coefficient 
f dischargefor the structure may be obtained with appropriate measurements 
f pressure, discharge, and valve or gate opening. This coefficient is common- 
y expressed as 


Q 
Cg ia iste aie ee ee tie aes 

d flat (5) 
1 which Q is discharge, H is head drop across the structure, and A is a de- 
ning area similar to that of the conduit, the control body, or the opening be- 
eath the gate leaf. Force measurements, such as downpull on a gate leaf, 
rque on turning vanes, or thrust on walls, may be obtained by applying mea- 
ured pressure heads on areas that will be affected, or by making direct force 
1leasurements. If the latter method is used, care must be taken to eliminate 
awanted model friction and extraneous vibration. 

The usefulness of air models in studying closed conduit systems has been 
iscussed. Such models discharge into the atmosphere and therefore operate 
submerged,” which is analogous to a water model discharging under water into 
large pool. Air-model data must therefore be interpretedas submerged flow. 
his is nohandicap in closed conduit systems, but must be recognized and ac- 
gunted for if model studies are made on prototype structure such as a gate 
at discharges water freely to the atmosphere.10 

Hydraulic models sometimes combine both open and closed conduit flow 
ynditions. An outlet works illustrates a case where studies are made of the 
mntrol valve or gate and of the stilling basin downstream. Both the open and 
losed conduit portions of the model are analyzed by the Froude relationships, 
the operating Reynolds numbers are in the appropriate range. 

It is difficult to build a model with equivalent dynamic response to that of 
e prototype. This does not preclude making studies of the stability of the 
ystem or reasonable predictions as to whether or not resonance and vibration 


{Ee Se ee SS 

8 “Cavitation Characteristics of Gate Valves and Globe Valves Used as Flow Reg- 
ations under Heads up to about 125 Feet,” by J. W. Ball, Transactions, ASME, Au- 
ist, 1957. 

9 “Cavitation Scale Effect,” by R. T. Knapp, Proceedings, Seventh Genl. Meeting, I. 
H.R., Lisbon, 1947, A6. 

10 “Air Model Studies of Hydraulic Downpull on Large Gates,” by W. P. Simmons, 
*,, Proceedings, ASCE, Vol. 85, No. HY 1, January, 1959. 
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will occur. Simultaneous recordings of pressure cell measurements at appro- 
priate points within a system can show if pulsations and/or surge will be pres- 
ent. These data may be applied mathematically to the prototype structure for 
analysis of possible resonance. In such tests, particular care must be taken 
to supply water to the model without pump or other mechanically induced pres- 
sure fluctuations. Pressure cell records will be requiredon the incoming pipe- 
line to show whether or not such pressure fluctuations exist. Special and de- 
tailed data analysis are required if such fluctuations are entering the system 
to complicate the pressure conditions. 

Graphs, charts, and other plots are usually needed to analyze and illustrate 
the data obtained. Final results find more usefulness when expressed in di- 
mensionless parameters that may be applied to other geometrically similar 
structures. Plots showing specific pressure conditions and other character- 
istics applicable to the structure under study are also made. 

Reports constitute the chief instrument through which the benefits of the 
model study are transmitted to designing engineers, and to the profession in 
general. They must be complete, accurate, timely, and to the point. However, 
brevity is no virtue if attained at the expense of completeness, clearness, OI 
accuracy. 


CONCLUSIONS 


1. An understanding of Reynolds similitude factors, and proper applicatior 
of them in designing, operating, and interpreting closed conduit hydraulic models 
will lead to reliable and accurate predictions of prototype performance. 

2. Models used in the usual studies of hydraulic controlelements and con- 
duit systems need not be excessively large, or operated with extreme flow ve- 
locities to yield reasonable model-prototype similitude. As a rule of thumb 
a model operating at or near a Reynolds number of about 1 x 106 is adequate. 

3. Surfaces of model flow passages must usually be smoother than the 
prototype surfaces for friction to be similar. In short segments or portions 0: 
conduit systems, full similitude in this respect is not important, and good re. 
sults will be obtained if the velocity distribution entering the section is turbu: 
lent and fully developed. In models with long water passages, more effort may 
be required to obtain friction similitude, if required, because friction would bi 
of relatively great importance. 

4. Particular care must be taken to produce accurately shaped flow sur. 
faces, and to maintain precise alinement of the model components. 

5. The selection, location, construction and maintenance of measuring da 
vices, such as piezometers and recording instruments, must be done with car 
and a knowledge of the type of data to be obtained. 

6. Tests should be continued to the point of obtaining all feasible informatid 
from each model so that the time and expense of the model construction will a 
exploited fully. 

7. Detailed diaries or logs should be kept of the studies, and comprehensiv 
reports prepared to make the material available to the engineering profession 
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HYDRAULIC MODELS OF THE ST. LAWRENCE POWER PROJECT2 
Closure by J. B. Bryce 


J. B. BRYCE.1—The author wishes to express his appreciation of the dis- 
ussions prepared by Mr. G. T. Berry and Messrs. G. B. Fenwick and F. R. 
3rown, both of which contained much additional useful information. As men- 
ioned previously, the comprehensive nature of the paper necessitated a gen- 
ral coverage of the subject and space did not permit treatment in detail of the 
arious problems encountered. However, in both discussions, interest was 
hown in the design of energy dissipation works and more information was re- 
uested in connection with the investigation of hydraulic loading of gates. 

With regard to the energy dissipation at the Iroquois Control Dam, Messrs. 
‘enwick and Brown note that only 10 of the 32 sluices were thus protected. As 
hey surmise, very substantial economies were effected by this partial treat- 
nent. A careful analysis supplemented by model tests indicated that only dur - 
ng construction would an energy dissipator be required, and in this high head 
riod the flow could be confined to 10 sluices. This was later corroborated in 
ctual operation. These discussers also note that the pattern of gate openings 
an have a considerable effect on energy dissipation and in the magnitude of 
he currents downstream. This was forcibly illustrated at Iroquois during con- 
tructionand even after the power pool was raised. Initially during these peri- 
ds, fairly severe currents were observed below the dam which proved trouble- 
ome to navigation. A very intensive series of model tests was run in the 
roquois Dam model, and a patternof gate openings was developed to maximize 
nergy dissipation and minimise downstream currents. When these patterns 
vere put in operationat the dam a remarkable improvement occurred and con- 
itions acceptable to navigation were produced. 

Messrs. Fenwick and Brown raised the question as to thenecessity of a 40- 
- wide end sill for the Long Sault Dam dissipator and whether a longer apron 
nd shorter sill might improve conditions. As mentioned in the paper, over 
nirty different designs were tested which included several combinations of 
1e above variables, and it was considered that the design selected was the best 
1 view of all the circumstances. Some of the factors which influenced the se- 
ction were: 1) the depth of the bucket was limited from geological consider - 
tions; 2) it was desired to have the apron as shortas possible, consistent with 
atisfactory performance; and 3) as the issuing velocities were still rather 
igh, the end sill was made wider than hydraulic conditions alone would require 
1 an effort to ensure that, if undercutting did take place, sufficient sill would 
emain to function properly. 

Both Mr. Berry and Messrs. Fenwick and Brown requested further details 


4 May, 1959, by J. B. Bryce. 
1 Hydr. Engr., Hydro-Elec. Power Comm, of Ontario, Toronto, Canada. 
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not permit a complete discussion of this subject, the writer feels that certain 
information additional to that given in the paper might prove useful. Accord- 
ingly, and with reference to the original paper, the following remarks are ap- 
pended. 

Port Gates at Long Sault Dam.—The results of the tests on these gates are 
given in Fig. 16 of the paper. The gates were 22 ft, 6 in. high, 17 ft wide, and 
2 ft 4 in. deep, and had the skinplate on the downstream side. The tunnel ports 
were 18 ft high at the upstream and 15 ft highat the downstream end and were 
16 ft wide. Each gate was to be capable of shutting off a maximum flow through 
the port of 14,000 cfs travelling at an average velocity of 40 fps and was to be 
fully submerged in operation. All tests on the gates werecarried out at a head 
of 77 ft on the sill and with free discharge conditions at the outlet. 

Two points not mentioned in the paper might prove of interest. The first is 
that when tests were carried out on the original gate design, the flat bottom 
caused the spring point of the flow to move from the downstream edge towards 
the upstream side under certain conditions. When this took place an oscillat- 
ing vibration of impressive proportions occurred. This behaviour is in agree- 
ment with that cited by Messrs. Fenwick and Brown in their discussion. The 
second point is the effect on the downpullof a recess made in the upstream face 
of the dam above eachport and immediately backof each gate when inthe wide 
open position. The recess in eachcase was 16 ft wide, 12 ft high and 1 ft deep, 
and it was found to reduce the downpull by 12 tons. It is believed that the re- 
duction in downpull by the recess is due to balancing the pressure across the 
top seal. However, if the recess was extended beyond the top of the gate in the 
fully open position, there was a tendency for the gate to move outwards from 
the dam. The significance of the clearance on the downstream side of the gate 
was also mentioned by Messrs. Fenwick and Brown. 

Following the tests reported in the paper, an additional investigation was 
made to determine the effect of several variables. In brief it was found that, 
without air venting, downpull varies directly with head, and that air venting is 
increasingly effective in reducing downpull as the head increases. It was also 
found that approach conditions have a considerable effect on the vertical forces 
in the gate. For example, lowering the bed upstream some 30 ft reduced the 
downpull by over 20%. 

Powerhouse Headgates.— The headgates tested were those for the Canadiar 
powerhouse which are of somewhat different design than those for the Unitec 
States powerhouse. Fig. 18 shows the gate dimensions and the details of one 
of the three water passages leading to the scroll case. The gate skinplate is 
on the downstream side and the gates operate in a fully submerged condition 
It was determined that the critical conditions for hydraulic loading on the gate 
would occur under runaway conditions with two of the inlets completely closet 
before the third gate started downward. The tests were performed with con- 
stant headwater and tailwater levels. The wicket gates were set to pass 14, 00t 
cfs at the beginning of the test (runaway discharge) and locked in position. 

Two of the passages were closed leaving the gate under test toclose the re- 
maining opening which was found tobe carrying about 13,000 cfs at 20 fps whe! 
the test commenced. In Fig. 18 the results show that the initial downpull of 1’ 
tons increased to a maximum of 119 tons at a gate opening of 8 ft. This may 
be compared toa gate weight of 37 tons. It was considered that the gate botton 
was sloped as muchas structural design would permit and no Poona loa 
was found possible. 
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Massena Intake Gates.—At the Massena Intake, two gates were tested, the 
amergency gate which had a skinplate on the downstream side, and the service 
sate which had the skinplate on the upstream side. Fig. 19 shows the arrange- 
nent and dimensions of the gates and the details of the water passages. The 
ests were carriedout at a constant headwater level and two selected tailwater 
evels. The maximum test flow was 21,000 cfs producing velocities of over 40 
ps through the port. 
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_ FIG, 18.—MEASURED VERTICAL FORCES ON CANADIAN HEADGATES 
FOR THE ST, LAWRENCE POWER PROJECT 


Emergency Gate Tests.—The emergency gate is upstream from the service 
ite and is completely submerged in operation. As shown in Fig. 19, maxi- 
um downpull forces of 68 tons and 76 tons were measured for the two heads 
sted. It may be noted that in the fully open position a downpull of about 50 
ns was measured, or about 140% of the gate weight. Comparing these downpull 
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forces with those measured on the Long Sault Dam port gates it would appear 
that they are somewhat greater than would be expected. It is believed that the 
two horizontal sections in the bottom, which accommodated the wheel, axles, 
increased the downpull force. Also, this gate is 50% higher than the Long Sault 
gate and the drag forces are probably greater. 

Service Gate Tests. —While the service gate is similar in size to the emer- 
gency gate, and the test conditions were identical, the results of the tests were 
greatly different as shown in Fig. 19. Downpull forces amounted only to 10 and 
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FIG, 19.—VERTICAL HYDRAULIC FORCES ON MASSENA INTAKE GATES 


15 tons respectively with the two heads, but uplift forces were encountered c 
16 and 10 tons respectively. A violent vibration was observed at and near 

gate opening of 25 ft at the higher tailwater level tested, when variations i 
force of plus or minus 35 tons were observed. This tendency was reduced bi 
not eliminated at the lower tailwater level. It is believed that the uplift wa 
due to the considerable back-pressure on the gate from the tailwater, and th: 
the vibration was due to unstable conditions at certain gate openings when 

rather violent pumping action took place in the tailwater. 
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Iroquois Sluice Gates, Long Sault Crest and Diversion Gates.—All these 
gates are similar in that they are 52 ft wide, 3 ft in thickness and have bottoms 
which sloped upward from the skinplate in the direction of flow at an angle of 
approximately 45°. All the gates had upstream skinplates and were unsubmerged 
in operation. The Iroquois gates were sheathed also on the downstream side 
io permit heating for winter operation. In Fig. 20 the results of all the tests 
and the gate dimensions and sketches of the gate settings are given. It should 
de noted that the head referred to is that on the sill for the Long Sault Crest 
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ate, and is the differential head in theother two cases. It was concluded from 
1e results that in the case of the Long Sault Crest gates and the Iroquois gates, 
e downpull was caused by drag forces only, as atmospheric pressure oc- 
irred underneath both gates at the point of maximum downpull. The increased 
ywnpull at the Long Sault Diversion gates was believed due to the restricted 
stting in which a great mass of water was thrown into the exposed downstream 
if 
2 
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side of the gate, thereby increasing the load. This was confirmed by testing 
this gate in the Iroquois setting where the downpull was greatly reduced. 

To obtain some information on the effects of bottom design on gates of this 
type the Iroquois gate was tested with a flat bottom. In the tests, at the point 
where submergence of the gate bottom began, a violent vibration occurred. 
The maximum downpull was also increased from 8 tons to 35 tons. 

While the tests previously described by nomeans cover the subject of ver- 
tical forces on gates, it would appear that certain general conclusions may be 
drawn. Vertical hydraulic forces have been shown to be very large at a num- 
ber of gate installations, particularly when the gate has a downstream skin- 
plate and is submerged. 

With a downstream skinplate, downpull is the dominant force which may 
many times exceed the gate weight. Substantial downpull forces may occur 
with the gate fully open if it is suspended close to a submerged port. Uplift 
appears to be relatively minor with this type of gate. The magnitude of the 
downpull appears to vary with the head on the gate and the size of the gate and 
is greatly affected by the shape of the gate bottom, being greatest with a flat 
bottom. While the major cause of downpull is undoubtedly the reduction in 
pressure under the gate due to the velocity, it has been shown that drag forces 
and approach conditions also have an appreciable effect. Effective air venting 
and equalizing pressures across a projecting top seal appear to reduce down- 
pull, but sloping the bottom members of the gate downward in the direction of 
the flow produces the greatest reduction. 

With upstream skinplates, downpull is not normally large, but uplift may be 
important if the gate is submerged. Sloping the gate bottom upwards in the 
direction of flow appears beneficial in reducing unstable operating conditions 
at part gate openings. 
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GROUND-WATER PROBLEMS IN NEW YORK AND NEW ENGLAND2 


Closure by Joseph E. Upson 


JOSEPH E. UPSON.1—The discussions of the writer’s paper by Messrs. C. 
3iemondand Robert O. Thomas reveal two quite different categories of ground- 
vater problems. Mr. Thomas’ comments elaborate on the need for adequate 
nvestigation of the geologic, hydrologic, and economic features pertinent to 
he development of ground-water resources. Solutions to these general prob- 
ems are certainly fundamental in the early stages of planning a water -supply 
nstallation, whether it be for irrigation, industry, or domestic supply. 

In contrast, Mr. Biemond points out a much more specific type of problem 
vhich lies in the narrower realm of the minute relationships between ground 
vater, with its dissolved, suspended, and bacterial loads, and the small-scale 
nvironment of its enclosing deposits. The effects of these relationships some- 
imes appear in operation of an installation and can be almost as critical as 
he amount of available water. Both categories are applicable and pertinent to 
levelopment of the ground-water resources of New York and New England. 
‘ew workers with, and students of, water are knowledgeable in the fields of 
oth categories. 

For the readers who may not fully appreciate Mr. Biemond’s finalreference 
9 advantages of artificial replenishment, the Amsterdam, The Netherlands, 
rater works has been pumping water from the Rhine within The Netherlands 
nd, with minor pretreatment, infiltrating it into a segment of the ground-water 
eservoir beneath The Netherlands coastal dunes since 1957. Apparently there 
as thus far been no clogging or other difficulties. 

The writer extends thanks and appreciation to Messrs. Biemondand Thomas 
wr interesting and informative discussions. 


———— anna EAE 
4 June, 1959, by Joseph E. Upson. 
1 Research Geologist, Ground Water Branch, Water Resources Div., U.S. Geol. Sur- 


y, Mineola, N. Y. 
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AN ENGINEERING APPRAISAL OF HYDROLOGIC DATA 


Task Group on Hydrologic Data of the Committee on Hydrology. 


DONALD W. VAN TUYL,! F. ASCE.—The thorough explanation by Mr. Kaz- 
mann of “estimated” and “safe” yields of ground water aquifers should greatly 
assist those considering the types of hydrologic data reviewed in the report of 
the Task Group on Hydrologic Data. Classifying “estimated yield” as interpre- 
tive data was intended to illustrate the distinction between quantities determined 
by compiling other values and quantities determined by more or less routine 
measurement or analysis. Primary emphasis was given throughout the report 
to basic and analyzed data with only bare mention of the interpretive classifi- 
cation. 

The twoterms which Mr. Kazmann discusses are usedin ground water com- 
outations without exact definition, as he clearly indicates. Both may be found 
io have specific values during engineering analyses, however, which may also 
invoive economic factors. 

The term “interpretive data” in the Task Group report should be viewed as 
an intermediate product of raw materials (data) but not an end in itself—a value 
or item of knowledge to be usedwith other items in obtaining a further answer. 
An item of datathat results from interpretation and extrapolation of basic data 
und is employed in subsequent hydrologic applications, therefore, may be 
classed as interpretive under this concept. 

For example, a computed flood stage for a particular frequency is the re- 
sult of numerous compilations of basic and analyzed data, but it is meaningful 
nly when applied as an item of data to an engineering problem. A member of 
he Task Group suggests that storage-yield curves for a surface water develop- 
nent are also interpretive data. 

One reference inthe Task Group report to another publication, unfortunately, 
yas out of date when the report reached printed form. Reference No. 1 should 
“ead: 


1. Langbein, W. B. and Hoyt, W. G., 1959, “Water Facts for the Nation’s 
Future,” The Ronald Press Co., New York; sponsored by the Conser- 
vation Foundation. 


? 


4 July, 1959, by Donald W. Van Tuyl. 
: 1 Chmn., Task Group on Hydrologic Data, Committee on Hydrology, and Water Re- 


ources Asst., Chamber of Commerce of U. S., Washington, D.C. 
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HYDRAULIC CHARACTERISTICS OF GATE SLOTS 
Discussion by C. T. Advani 


C. T. ADVANI.1—The writer found Mr. Ball’s paper most interesting, es- 
pecially concerning the large number of model test results that he has been 
able to furnish on various different type of slots. 

A certain number of the results with the different type of slots obtained by 
Mr. Ball agree with those found during the detailed program of investigations 
carried out during the last 3 yr in the laboratories of the Societe Grenobloise 
d’Etudes et d’Applications Hydrauliques, Grenoble. The investigations were 
especially concerned with the problem of partial openings of vertical lift gates 
under very high heads and the resulting cavitation conditions in the gate slots. 

The writer would, however, like to stress that it is absolutely essential that, 
in order to interpret the test results obtained correctly, the similitude between 
the prototype and the model must be established perfectly. 

If this is not done, in a number of cases the interpretations may prove to 
be rather faulty. Very great care on this account has to be taken in designing 
the test installation, so that the Froude similitude and the Reynolds number are 
sufficiently high. The test installation used in the laboratory had the clear 
gate dimensions: 600 mm high and 260 mm wide. The gate model was hy- 
draulically operated by servomotor control. The cavitation circuit was de- 
signed to pass 1200 liters of water per sec. With this type of circuit having 
satisfactory pressure control, it was possible to obtain reasonably accurately 
the approach of cavitation breakdown. Also, in order that the results obtained 
are not faulty, it is essential that the circuit be clear of parasite cavitation in 
bends, the pump, etc., so that the observed changes in the regime are due to the 
gate operation and hydraulic conditions only. Further, if water with some air - 
content is to be used, the circuit should be provided with air reabsorber, so that 
bubbles do not collect in dead pockets. 

The writer completely agrees with Mr. Ball that the slot problemis essen- 
tially a detailed study of pressure intensity and distribution. The heads, the 
dimensions as well as the downstream hydraulic conditions will determine 
whether there is cavitation or not. The intensity of cavitation will at the same 
time depend on the gate opening. For different hydraulic and slot conditions, 
the critical openings vary. For a number of existing installations, the gates 
have not shown signs of cavitation, because they have never functioned under 
the critical openings for the hydraulic conditions imposed. 

_ The writer feels that the correct mode of attacking the problem was not to 
find a typical design which would be fool-proof under all different operating and 
hydraulic conditions but to study each particular case separately. 


@ October, 1959, by J. W. Ball. 
1 Engr., Sogreah, Grenoble, France. 
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The series of tests that have been carried out by the writer on different 
types of slots and different hydraulic conditions imposed have proved that for 
effective heads of over 40 m with downstream undrowned (free discharge) and 
with typical classical vertical lift gates, freedom from cavitation cannot be had 
under all possible openings. This condition is also applicable for arrange- 
ments having upstream contraction and air duct tubes feeding along the height 
of the slots. Actually, the best arrangement is with upstream sealing and nor- 
mal slots without offset. With this disposition, it would be possible to operate 
the gates without large cavitation risks for higher heads, but only for inter- 
mediate (between 1/3 and 2/3) and full gate openings. 

If however the downstream jet is drowned (backpressure), the results are 
most unfavorable whatever the position of the sealing, the type of the slot, or 
the gate opening. 

If the vertical lift gates are required to function partially open under very 
high heads, then a solution other than the classical gate has to be envisaged. 
During our investigations, it was specified that a solution should be foundwhich 
should not only be suitable from the hydraulic point of view, but had to be feas- 
ible mechanically and from the economic stand point, the cost had not to ex- 
ceed those of the classical vertical lift gate. The solution of baffles or exten- 
sions to bottom corner of gate leaves (in any case not a very healthy solution) 
or very large offsets in the slots had to be rejected. After a considerable 
amount of study and tests, a gate design was evolved which satisfied the differ - 
ent conditions imposed. 

A brief mention of this solution was made by Messrs. Zaky and Sauvage de 
Saint Marc at the last Congress of I.A.H.R. 

It is impossible to give a complete idea of this solution in this short note of 
comments. Briefly speaking, the gate leaf knife edge is designed to be up- 
stream of the gate slot, it being supported on a sort of parabolic bridge struc- 
ture. The slots from the hydraulic point of view are nearly in uniform flow 
without large variations in pressure along the height of the slot. The results 
that were obtained during the tests were most encouraging. 

In case of partial gate openings, just downstream of the gate, one has over 
the resulting jet surface either air (free surface flow) or a mass of water 
(drowned jet) in which the pressure distribution is practically hydrostatic. In 
case of free surface downstream, the h downstream (back pressure) is the 
height of the jet itself above the gate sill. In case of drowned jet, the h down- 
stream is the height of water in the air duct, also based on the sill level of the 
gate. 

To obtain accurate results, very delicate pressure measurements have tc 
be taken. An exact idea of pressure fluctuations is essential. To know the 
average pressures is not enough. 

The cavitation coefficient can be described as: 


_ h downstream + h atmospheric - hy 
h upstream - d downstream 


The numerator is the difference between the absolute pressure downstream o: 
the gate and the vapour pressure hy. The denominator is = of the jet, whicl 
is in itself the difference between the pressure just on the upstream and down. 
stream side of the gate. The sill is taken as a datum for all heads. : . 
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The classical vertical lift gates with upstream sealing and normal slots 
fithout offset and downstream undrowned cavitate around o = 0.25 at difficult 
penings (less than 0.2 and greater than 0.7) and at o = 0.1 at favorable open- 
ngs. With downstream jet drowned, the o is rarely under 1. 

The writer trusts that the comments give at least a rough idea of the way 
he problem of gate slots has been studied and that they will help to compare 
yith the results obtained by Mr. Ball. 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL CHANNELS4 


Discussions by C. Blanchet, D. R. Dawdy and R. W. Carter, M. Gamal 
Mostafa, and James K. Culbertson and Carl F. Nordin, Jr. 


C. BLANCHET.1!—This is an extremely interesting investigation, from the 
ractical point of view, because it provides engineers with an easy homogene- 
us calculation method that covers nearly all the experimental results known 
sday, which is not the case with a certain number of earlier investigations. 
ike earlier publications, this article shows that the problem is not only very 
omplex, but also that despite the many experimental results that are known 
nd applied, their number and organisation are still grossly insufficient. 

In spite of this well known weakness, that affects us all, the authors—like 
1any of their predecessors—have given in to the temptation to generalise the 
olution and to give it a vaguely theoretical appearance. They have, for in- 
tance, considered parameters for which little experimental data, or none at all, 
reavailable. The results are expressed in terms of such parameters as grain 
hape, grain and fluid specific gravities, viscosity, gravity and temperature, 
one of which have really been considered in the tests under analysis. Some 
oubts therefore remain as to the reliability of the results, if these are pre- 
isely the parameters to be applied. 

For instance, although the introduction of the drag coefficient Cg deduced 
rom the fall velocity of a grain in still water is probably valid for the more 
r less spherical grains used in the tests, it is doubtful whether it would still 
pply, if, say, pins were used instead. Under these conditions, the form effect 
fan isolated pinunder free fall conditions would probably bear no relationship 
» that of a pin lying flat on the surface of the water. 

Considering a sphere, it is also doubtful whether the sudden change in the 


rag coefficient occurring when 


nder free fall conditions incalm water is also valid for a bed of spheres. The 
ifference betweenthe state of the flows and the geometry of the system is too 
reat. In the writer’s view, the introduction of the fall velocity of an isolated 
rain is not representative of the form effect of the same grains when forming 
bed. 

A further doubtful result, which is probably due to the fact that the density 
f the grain has not been considered over a wide range, is that the apparent 
eight of the grain in water does not affect the velocity under turbulent con- 
itions (Eq. 26). In other words, for the same rate of flow, the same slope, the 


2 November, 1959, by Hsin-Kuan Liu and Shoi-Yean Hwang. _ 
1 Chf. Engr., Sogreah, Grenoble, France. 
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same canal, and grains of the same diameter but of widely varying specific 
gravities, (with a high sediment load in one case and none in the other) the ve- 
locity is supposed to be the same. Is this really so? Or could it be that this 
result was obtained because the tests considered were not carried out over suf. 
ficiently wide specific gravity and sediment load ranges ? 

As the title of the article suggests, it is considered that the results of the 
investigation are very true for natural alluvial channels, in which such parame: 
ters as gravity, grain and liquid specific gravities, grain shape, viscosity an 
temperature only vary slightly. However, contrary to the apparent possibili. 
ties shown in the presentation of the results, some doubts arise as soonas ont 
leaves the scope of the basic experiments. According to the author’s conclu: 
sions, it should be possible to settle the question once and for all by carryin; 
out tests on a very wide scale. 

The writer feels, however, that such things as grain shape, specific gravity 
sediment load and viscosity should be investigated very thoroughly. Testi 
should be carried out with a variety of shapes like pins, discs, lenticular forms 
and so forth, with relative densities varying between 1 to 100, and with ver: 
high sediment loads. 

Inthis connection, we would like to draw attention to anarticle by M. Rottne: 
in “La Houille Blanche,” n° 3 (June, 1959). M. Rottner collected 2,500 experi 
mental results and succeeded in classifying them in terms of the followin 
parameters: 


in which G is the sediment load, and the remaining symbols are as in the arti 
cle under discussion. 

M. Rottner’s investigations not only enables the same determination to b 
made, but also gives the sediment load. It is thus more complete and opens u 


FIG, 1.—DATA FOR STRAIGHT ALLUVIAL CHANNELS 
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new prospects on a certain number of peculiarities and on zones so far inade- 
quately investigated by empirical methods. The graph obtained, Fig. 1, shows 
certain peculiarities that will readily be noted by specialists in the subject. 


DER. DAWDY,2 and R. W. CARTER, ? F. ASCE.—Mr. Liu’s discharge for- 
mula for alluvial channels is ingenious, but his basic assumption that various 
empirical flowformulae can be synthesized into onerational form does not ap- 
pear to be justified. This is one universe to which order must be brought not 
through synthesis but through catharsis. 

The author used the Blasius formula for grain sizes smaller than 0.1 mm. 
The Blasius formula is valid only in the range of Reynolds numbers less than 
80,000, which means that for v = 1.12 x 10-5 for 65° (as used by Liu) the fol- 
lowing limitations hold, where R = 4 V D/V 


R = 4,000 R = 80,000 
V, in D min, D max, 
ft per sec in ft in ft 
10 -0011 .022 
5 -0022 .045 
2.0 .0045 .09 


This range of variables is not one often encountered in nature for the flat 
smooth bed regime. In fact, it probably is not encountered often even in labo- 
ratory flumes. In nature the flat smooth bed regime seldom occurs at depths 
less than half a foot, and then only for velocities on the order of 5 fpsor great- 
er. Thus the Rmin should be on the order of 8 x 105. Actually from field ex- 
perience, the smallest R known to the writers is on the order of 18 x 109. 

Perhaps the major appeal of the Blasius formula is thatit is a simple power 
function, and easily handled mathematically. Thus, simplicity seems to be the 
mother of invention. Other empirical formulas? have been offered as well as 
the Blasius to describe the f @ R relation just beyond the critical R. Blasius 
considered a sufficiently small range of R so that a first order approximation 
could be used. Since in actuality the relation is not a simple power function, 
the range of R can be increased only by using a more complex equation. Her- 
man proposed, 


This equation would be more difficult to use than the Blasius and is inadequate 
in that it covers only the lower limit of field data. Thus, either formula is at 
best only an approximation, and perhaps ‘should not be used for determining the 
powers of R and S for flow in open channels. 

The discharge formula proposed by Mr. Liu was used to compute the ve- 
locity for three reaches of stream channel where measurements of velocity, 
energy slope, and depths of flow were available. A comparisonof the observed 
and computed velocities in the dune and smooth flat bed regimes is shown on 
Fig. 1(a). The measured velocity for all measurements is higher in the dune 
regime and lower in the flat smooth bed regime than the velocity predicted by 
the Liu equation. The range of scatter is on the order of 500%. 

2 Hydr. Engr., U. S. Geol. Survey, Washington, D. C. 


4 


3 Hydr. Engr., U. S. Geol. Survey, Washington, D. C. J 
_ 4 «Rssentials of Fluid Dynamics,” by L. Prandtl, Hafner Pub. Co., New York, 1952. 
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A similar plot for three additional reaches is shown on Fig. 1(b). For these 
aches, the computed velocitiesfor all flat smooth bedregime measurements 
e too high, and those for all but one dune regime measurement are too low. 
ne Slopes usedin computing these velocities from the Liu equation were taken 
om topographic maps. The predicted velocity appears to be a better approxi- 
ation of the measured velocity for the larger grain sizes. 

The criterion offered for prediction of regime (Fig. 1 of paper) is based 
yon V,/w andwd/v. This plot first was introduced ina closure to a paper 
‘Mr. Liu and therefore has not received full discussion. Wd/v has replaced 
,a/v in a plot Liu previously offered. Mr. Liu states, “This is permissible 


| eS a a 


oElkhom, dune regime 
xElkhorn, plane bed regime 


e Bernalillo dune regime 
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FIG. 2.—COMPARISON OF MEASUREMENTS OF KNOWN BED CONFIGURATION 
WITH LIU’S PREDICTED BED CONFIGURATION 


cause V, d/y = S ; woe Actually this is permissible only within the mean- 
‘of dimensional analysis. Even there it is not entirely permissible, since w 
idare not independent. “The mechanics upon which the criteria for vari- 
3 bed configurations are based . . .” does not®enter into either form except a 
steriori. Actually, there are six possible Reynolds numbers combining V, 
,and w with R and d. The four involving V, and w are identically acceptable 
pirically and equally unacceptable if a physical significance is implied. The 
m w d/y is chosen for ease of use only. This field is a muddy one, and 


5 “Mechanics of Sediment Ripple Formation,” by H.K. Liu, Proceedings, ASCE, Vol. 
I HY 2, April, 1957. : 
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mere manipulations of dimensional analysis will not make it any clearer. / 
any rate, the criterion offered using v,/w and w d/y does not define the be 
configuration of natural channels. Data for the Rio Grande River near Be! 
nalillo, New Mexico, and for Elkhorn River near Waterloo, Nebraska, are shov 
in Fig. 2 to emphasize this point. 

Mr. Liu further states that the Froude number “can be considered as 
index indicating the energy loss caused by the surface waves.” Since 


the implicationis that for a given measuring site in the field where the ener; 
slope remains approximately constant, the energy loss caused by surface wave 
is a function of the Chezy C. For the flat smooth bed regime, Chezy C is ¢ 
the order of twotimes that for the dune regime, and is essentially equal to th 
for the antidune regime. Therefore, it follows that the energy loss caused | 
the surface waves is about the same for flat smooth bedand antidune regime 
The plane bed has no surface waves, the antidune violent surface waves. Als 
it is implied that the energy loss due to waves for the flat smooth bed regin 
is twice that for dune regime, although for the flat smooth bed regime bo 
water surface and bed are plane. It is believed the physical significance a: 
signed to the Froude number may be more apparent than real. 

The author should state whether the points used to define Cg in Fig. 11 a: 
individual points, or are averages of a series of runs. From the small amou 
of knowledge gained through working with datafor alluvial streams, it does n 
seem that the data ever plot as well as indicated on the figure. The scatter | 
Figs. 1 and 2 herein seem more typical. 


M. GAMAL MOSTAFA.®—The authors have tackled an important proble 
in hydraulic engineering and have presented a thorough study of the prese 
status of research towards its solution. Attempts to extend the applicabili 
of flow equations, which were originally derived for rigid boundaries, tot 
field of alluvial streams have so far been unsuccessful. The authors ha 
therefore turned to empiricism and have proposed a flow formula, Eq. 31, co 
taining three sediment factors x, y and cg which vary for different bed mate 
al and configuration. As a guide in the application of their method, they ga 
Fig. 1 which can be used when needed to predict bed configuration. 

To estimate x, y and c, the authors gave Figs. 9, 10 and 12, which were d 
duced from some collected data. The considerable jumps in the values of ¢ 
x and y when sand ripples appear or dunes disappear may give the impressi 
that suchchange in bed configuration occur suddenly. If this was true, it 
thought that there should be a critical state for such a sudden transfer. Co 
sidering the movement of a sand bed under an increasing flow traction, | 
change of bed shape from plane to ripple is usually noticed just after the k 
ginning of general bed movement. The development of sand ripples is gradi 
and their transformation into dunes is also gradual. This was first illustra’ 
seventy years ago, by Deacon, 7 who also noticed the gradual flattening of | 


6 Dir., Hydr. Research Sta., Delta Barrage, Egypt, V.A.R. ~ a 
7 “An Introduction to Fluvial Hydraulics,” by S. Leliavsky, Constable and Co., L 
London, 1954, p. 11. 
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ed at relatively high flow velocities. L.G. Straub,8 found that the maximum 
eight of dunes was reachedat critical flow conditions, corresponding to mini- 
lum specific energy, and that as the velocity of flow rose above the critical, 
¢ undulations became less and less until finally the bed became smooth. 
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‘ FIG, A.—NILE WATER PROFILES BETWEEN ASWAN AND ISNA 

“Fig. A gives the water surface profiles of the Nile between the Aswan dam 
id Isna barrage at six different discharges ranging from very high to very 
Ww. During the field measurements, in each case, the discharge was kept 


8 “Discussion on Sand Movement in Fluvial Models,” by L. G, Straub, Transactions 
+E, Vol. 61, 1935, p. 867. 
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constant for a period long enough to maintain steady flow. A certain more 0 
less straight reach, 30 km long, located upstream of the backwater effect fe 
Isna barrage, was chosen for testing the authors’ flow formula. Mechanica 
analysis of bed material in that reach gave a mean sand diameter of 0.25 mm 
Table 1 is a summary of the results. In the computations x and y, which wer 
deduced from Figs. 9 and 10 (dune bed), were taken as 0.44 and 0.31 respec 
tively. The flood discharges corresponded to the transition region in Fig. 
and their computed cg were 11.61 and 11.40. For medium discharges cg wa 
reduced to 10.1 and the bed configuration was “dune.” At low discharges c 
rose again and reached a value of 16 at the lowest discharge. This was ex 
plained by the authors’ Fig. 1 which showed that the bed condition was close 
to “ripple.” 


TABLE 1.—APPLICATION TO NILE DATA 


<a ae Ses Mil lea ee a ee TE ee ] 7 a 
Discharge] Discharge water slope | V, = VegDS 

Q, incum|] per unit depth 
per sec width q, D, in x 10° |cms per sec 

in mper sec | metres 


6 
(1) (2) ; (3) (4) (5) ( 4 
2.6 


10900 15.57 9.5 7.4 8.3 2 
10200 14.57 9.3 7.0 7.95 2.5 
7500 10.71 8.33 6.3 7.16 2.24 
3360 4.80 4,90 5.8 5.27 1.64 
1620 ona B72 5.6 3.87 sie? 

STL pe “hued 1.42 | 5,1 2.66 .83 


Application of the authors’ formula to the Nile data has therefore prove 
satisfactory. Discharges estimated by the formula differed from measure 
discharges by no more than 20%. Bearing in mind that Manning’s formula, fo 
example, may lead to an error of over 100% in the estimation of river flov 
the authors’ formula can certainly be considered a great improvement. 

It would be very interesting to know how the proposed formula fits in it 
application to the Missouri River, which is knownto carry a heavy load in sus 
pension, and also in its application to alluvial streams with steep slopes. 


JAMES K. CULBERTSON,9M. ASCE, and CARL F. NORDIN, Jr.10_The ar 
thors are to be commended for an analysis which covers a wide scope of tt 
available flume and canal data, and for the presentation of a simplified formu 
which, if proven applicable, may be adapted withease to either design problem 
or to the prediction of channel behavior. Moreover, the relationship betwes« 
bed configuration and mean velocity of flow appears to be related to the phe 
nomena of stage-discharge discontinuities which exhibit themselves in true a 
luvial channels, and the results of the application of the proposed formula 
data from measurements of the Rio Grande in New Mexico suggest that tl 
formula may be adapted to some natural channel conditions. However, certa 
modifications are necessary due to the fact that the authors’ Fig. 1 does n 
adequately define bed configuration in the field case. 

Current studies (1960) indicate that the true sand channel reaches of the R 
Grande exhibit stage-discharge discontinuities which apparently represe 

9 Hydr. Engr., U. S. Geol. Survey, Albuquerque, N. Mex. 

10 Hydr. Engr., U. S. Geol. Survey, Albuquerque, N. Mex. 
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anges in regimes of flow, and which may be detected by plots of hydraulic 
dius versus velocity. (A true alluvial or sand bed channel as used in this 
scussion has the same characteristics as the authors’ alluvial channel, in that 
¢ bed material is granular and cohesionless and the amount of supply of sedi- 
ent is equal tothe amount of sediment transport.) Fig. A represents one such 
ot of selected measurements of the Rio Grande near Bernalillo, N. Mex. The 
yw may be divided into a lower and an upper regime, with a transition regime 


Hydraulic Radius in Feet 


Velocity in fps 


FIG. A.-RELATION OF HYDRAULIC RADIUS TO VELOCITY FOR RIO 
GRANDE NEAR BERNALILLO, NEW MEXICO 


ween. The exact nature of all changes in bed configuration which occur with 
inges in regimes of flow is not known, but because of the decreases resist- 
e toflow which must accompany these changes, it is assumed that the lower 
rime represents flow over dunes and sand bars, and that during the tran- 
. to the upper regime, these bed features are obliterated or flattened by 
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the actionof the flow. It was found that the velocities for points which fall alor 
the lower regime curve may be computed using the authors’ curves for the sar 
bar configuration, and that velocities for points in the transition stage betwee 
lower and upper regime may be computed using the proposed curves for tl 
flat bed configuration. Velocity measurements which plot along the upper re 
gime curve in Fig. A appear to represent flow conditions associated with a fl: 
bed or antidunes and standing waves, and for this regime of flow it was poss: 
ble to modify the Cg curve to actual field conditions, as suggested by the at 
thors. 

A total of 73 measurements was applied to the proposed formula. The: 
measurements covered a range of depthfrom about 1 ft to9 ft, of velocity fro 
2 fps to 8 fps, and of median particle size of bed material from 0.14 mm to 0. 


mm. Values of we and v,/w were computed using average median partic 


size and average water surface slope for the reach. These values were plo 
ted on the authors’ Fig. 1. Eleven points fell in the dune bed configuration a1 
26 points fell in the transition (flat) bed configuration. Of the 37 points, only 
of the values computed using the authors’ Ca, x, and y were within 20% erro 
Classification into regime then proceeded on the basis of plots of hydraul 
radius versus velocity and of field observations made at the time of the me: 
surement of water surface appearance andof relative bed compaction. Veloc 
ties for these points were then computed using the authors’ values of Cg, x, al 
y, and assuming that the lower regime points represent dune bed configurati 
and that the transition regime points represent the flat bed configuration. Abo 
60% of the values so computed were within 20% error, showing a marked in 
provement over the classification of bed configuration based upon v,/w a 
wd 

pale 

In general, the field observations were of value in indicating where stan 
ing waves were present and where the channel was hard and flat, which in < 
cases were associated with the upper regime. The break in the hydraul 
radius-velocity curves was not so well defined for all other stations as for t 
Bernalillo station, and considerable difficulty was encountered in differentiati 
the lower and the transition regime, especially where observations of relati 
bed compaction and water surface appearance were not noted. 

For all other measurements which fell in the upper regime, a value of | 
was computed using the given exponents for x and y for a plane bed and t 
measured velocities, and average values were used to define the modified | 
curve shown in Fig. B. The lower portion of the curve falls considerably b 
low the values which might be expected for these particle sizes, and it was d 
termined that for these stations portions of the bed scoured down to rock a 
gravel at higher stages of flow. Further, the writers believe that reaches 
the Rio Grande which have median particle sizes of bed material greater th 
about 0.40 mm cannot be classifiedas true sand channel streams; that is to s: 
the bed configuration, velocity, and roughness parameters show no unique 1 
lation with change in regime of flow. This fact limits the usefulness of | 
curve to those reaches which meet these criteriafor true sand channel strean 

The results of this study indicate that the proposed method for comput 
mean velocity of flow can probably be applied successfully to true sand chz 
nels such as portions of the Rio Grande, but that additional criteria beyond | 


d 
relation of = to V,/w are needed to define adequately bed configurations. 


) as Pla rae 


fi 


DISCUSSION 


SHMRBANGE 


TTA 


juaoyjeod ebsoyosid 


30) 


rans pan 


a 


mm 


in 


~ Median Size d 


FIG. B.—VARIATION OF C, WITH d FOR FLOW OVER ALLUVIAL BED 


102 June, 1960 HY 


The writers further believe that any criteriafor determining regime of flo 
in true sand bed channels that is based on shear velocity cannot be sufficien 
If a constant energy gradient is assumed for a reach which has the hyaraul; 
radius-velocity characteristics indicated in Fig. A, then shear velocity be 
comes a function only of depth of flow or hydraulic radius, and it is possib. 
that more than one flow regime can exist for a given hydraulic radius. 

In conclusion, the writers propose that the following items be included | 
the authors’ suggestions for future research: 


1. The method should be extended to include the “antidune” bed configuratic 
as additional information becomes available. 

2. Additional field data are needed, with special attention directed to tt 
development of equipment and techniques for measuring actual bed configt 
rations. River profile studies such as the one conducted by W. C. Carey a1 
M. D. Keller11 shed considerable light upon the changes in bed configuratior 
which occur with change in stage, and further investigations along this lit 
would be of great value. 

3. The relation of sediment discharge to bed configuration, and its effe 
upon the proposed formula should be investigated. 


11 «systematic Changes in the Beds of Alluvial Rivers,” by Walter C. Carey and 
Dean Keller, Proceedings, ASCE, Vol. 83, No. HY 4, August, 1957. 
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EFFECT OF AQUIFER TURBULENCE ON WELL DRAWDOWN2 
Discussion by E. G. Kruse 


E. G. KRUSE.1-—The effect of turbulence on head losses in aquifers near 
umped wells has been considered by the author. In the West many wells are 
onstructed with an envelope of gravel to prevent sand from moving into the 
ell. In the light of this paper it would be interesting to consider the effects 
f turbulence on flow in gravel packs. 

In the gravel pack both velocity and grain size are larger than in the aqui- 
1, therefore turbulent flow is more likely to exist. The permeability of the 
ravel pack should be several times higher than that of the aquifer. The head 
SS resulting from turbulence, therefore, may not be important. 

A series of tests have been made by the writer for the purpose of determin- 
ig proper size and uniformity for gravel packs for use with wells in various 
juifers. The test model represented a 30-in. diameter well with a 12-in. 
ameter screen, a size representative of many irrigation wells. Seven pie- 
ymeters along a radius in the gravel pack permitted head readings. Dis- 
larges ranged from the equivalent of 10 to 100 gal per min per vertical ft of 
ell. Reynolds numbers were computed by the method used by the author. 
eynolds numbers in the gravel packs at the well screen varied from 2.3 to 
1.0. Most were greater than 10. 

To determine whether flow in the gravel pack of the model was in a lami- 
ir or transition regime, the Reynolds number at the well screen was calcu- 


ted. The value of n in the expression V" = K 4 was determined according 


‘the author’s table. Then using the head readings at two points in the gravel 
ick, one near the pack-aquifer interface and the other near the well screen, 
evalue of K was calculated for two cases, n as listedunder the heading “Con- 
usions from Analysis of Test Data” and n equal to 1 (assumed laminar flow). 
he head loss curve in the gravel pack was plotted for each assumed n and 
mpared with the experimental data for intermediate points in the gravel pack. 
For Reynolds numbers on the order of 20, the assumed laminar flow curve 
rrees slightly better with the experimental data than does the transition curve. 
1is is illustrated in Fig. 1 for one test typical of several at this Reynolds 
mber. The agreement of the data with the laminar curve is therefore ap- 
rently real rather than due to experimental error. For tests with higher 
3ynolds numbers at the well screen the experimental data generally plotted 
tween the assumed laminar and transition curves. It may be concluded that 
me head loss due to turbulence often exists in gravel packs for Reynolds 
: 
2 November, 1959, by Joe L. Mogg. 
1 Agric. Engr., Western Soil and Water Met. Research Branch, Agric. Research Ser- 
ce, U. S. Dept. of Agric., Fort Collins, Colo. ; 
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numbers greater than 20. The values of velocity exponent listed by the autho 
are probably slightly high due to the factor of safety that the author has includec 

The head losses in the model tests were 6 in. or less through properly se 
lected and developed gravel packs. This value would include the component c 
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FIG, 1.—HEAD LOSS CURVES FOR ASSUMED LAMINAR AND TRANSITION FLOWS 
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FIG, 2.—HEAD LOSS THROUGH WELL SCREENS 


loss due to turbulence. For a pumping lift of 50 ft or more the total head los 
through the gravel pack would amount to less than 1% of the total lift and tt 
component due to turbulence could surely be neglected. 

Reynolds numbers in the aquifer at the pack-aquifer interface were most! 
less thanten and inall but two tests were less than 20. It would seem that he 
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oss due to turbulence in the aquifer of a large diameter gravel-packed well 
‘an also be neglected in computing drawdown in wells. 

Experimental data obtained byG. L. Corey, Jr.2 (Fig. 2) support the author’s 
bservation that screens with insufficient open area may cause significant head 
osses. Corey’s tests showed that head losses through the screen decreased 
iS Cc, the percentage of open area, increased. For Cg above some critical 
ralue, head losses remained constant. This critical value was dependent on 
lischarge and the size of material surrounding the screen. Corey also found 
hat the head losses through the gravel packs were quite small compared to 
otal losses through the well screens. 


a eS a ee ae 
2 “Hydraulic Properties of Well Screens,” by Gilbert Lee Corey, Jr., Unpublished 
faster’s Thesis, Colorado State Univ., June, 1949. ; 
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THE VORTEX CHAMBER AS AN AUTOMATIC FLOW-CONTROL DEVICE2 
Discussion by John C. Stevens 


JOHN C. STEVENS,! F. ASCE.—This paper was obviously submitted before 
‘aper 1461 of the author’s reference 5 was printed in Transactions as Paper 
004, Vol. 124, 1959, p. 871. In Transactions, the Vortex Number of the au- 
iors’ signs and symbols has been changed tothe Kolf Vortex Number and given 
1e symbol K in recognition of the original contribution to the science of hy- 
raulics by Richard C. Kolf, M. ASCE. In this discussion the symbol K is, 
ierefore, used instead of V. 

The authors have added to our knowledge regarding vortex flow through 
orizontal orifices. In particular they have shown that the orifice discharge 
oefficient C can be determined approximately from the velocity of the water 
ntering the vortex chamber. They have also established the relation between 
and @, the angle of divergence of the water from the horizontal as it leaves 


1e orifice, to be 
A KO. COS Wj ened ae pheeieeis + nea (1) 


therwise K had heretofore been determined from the circulation I, the orifice 
iameter D and the average velocity of water, through the orifice, to be 


epee GRC Ce TIE MARES SEBLN -h (2) 
= D/2gH 


As stated by the authors, the water surface profile computations were “based 
n a measurement made only one diameter away from the orifice” that is, by 
ubstituting D for X in Fig. 1. It was found in the University of Wisconsin ex- 
eriments that coordinates of this point were a fair index of the hyperbolic 
ater surface profile of a free vortex? given by (28) except perhaps near the 
rifice where viscous forces are in evidence. 

The writer was interested in knowing how closely the values of K by the two 
ethods agreed. From the authors’ table has been prepared. 


K/ cos 


.89 
son 
1.08 
97 


2 December, 1959, by R. C. Kolf and P. B. Zielinski. 

1 Cons. Hydr. Engr., Stevens and Thompson, and Leupold and Stevens Instruments, 
C., Portland, Ore. 

2 Transactions, ASCE, Vol. 124, 1959. 
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The average of the ratios in the last columnis 0.985 indicating that the aver- 
age value of K by either approach was within 1-1/2% of each other, although 
the maximum variation was 11%. With a more refined method of measuring 
this angle, which the authors admit was difficult in the Marquette University 
laboratory, better individual agreements could doubtless be expected. 

A comparison of the discharge coefficients versus the Kolf vortex number 
for both the Marquette and the Wisconsin experiments are shown in Fig. 1. 
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FIG, 1.—COMPARISON OR COEFFICIENTS AND VORTEX NUMBER 


Curve M is takenfrom the authors’ Fig. 2 and is defined by experimental datz 
while Curve W is from the composite curve of the authors’ reference 5, th 
solid portion of which is similarly defined. The deflection toward Grandma’ 
near the bottom of Curve M might disappear if some allowance could be in 
cluded for the non-hyperbolic portion of the vortex surface profile as wate 
approaches the orifice. 
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THE SETTLING PROPERTIES OF SUSPENSIONS? 
Diseussions by Lucien M. Brush, Jr. and Hau-Wong Ho, and E. J. Hall 


LUCIEN M. BRUSH, Jr.1 and HAU-WONG HO.2—The author deserves con- 
siderable praise for his lucid presentation of an interesting subject. The de- 
ailed presentation of the continuity equation for a sediment-water suspension 
inturbulent motionis straightforward and pertinent to the ensuing experiments. 
The experiment and the interpretations of the results are clearly presented, 
as are the comments on similarity, and the speculations as to the application 
of the multiple-depth analysis to sediment problems involving open-channel 
‘low and sediment removal. In addition to the interpretations made by the au- 
thor concerning the experiment, several other interpretations exist which the 
writers feel are worthy of comment. Because the core of the author’s paper 
is based on the so-called measurement of w, the local mean settling velocity, 
this discussion will be directed toward the experimental and interpretive as- 
sects of the investigation. 

In presenting the continuity equation for sediment suspensions, the author 
states that turbulence will appear explicitly through the diffusion terms, where- 
4s flocculation or hindered settling due to turbulence will occur implicitly in 
the fall-velocity term. Nevertheless, the results obtained by the author for 
the multiple-depth analysis are interpreted to be due to flocculation alone. 
Thus, any effect due to turbulence, either implicitly or explicitly, is discount- 
2d. There seems to be little justification in ignoring turbulence, because it is 
most certainly present, at least initially, if the procedure outlined by the author 
is followed in making the analysis. Turbulence enters the system as a result 
of “thoroughly mixing” the sample (to obtain a uniform dispersion) and as a re- 
sult of pouring the sample into the measuring cylinder. Obviously, the intensity 
of the disturbance caused by this experimental procedure will decrease with 
ime due to the viscous effects of the fluid, but the net outcome of the intro- 
juction of decaying turbulence will be to generate curves similar to those ob- 
ained by the author in Fig. 8—with or without the presence of flocculation. 

Although it is difficult to measure the scale and intensity of the introduced 
urbulence without running rather elaborate experiments, it is possible to make 
1 qualitative estimate of the amountof time that the turbulence effect would be 
of importance in an experiment similar to that made by the author. The writers 
jlaced fluid droplets with a specific gravity slightly greater than one in a con- 
ainer filledwith water and thoroughly mixed the liquid to obtain a dispersed 
suspension. The contents of this container were then quickly poured into a 
sylinder which had dimensions similar to the one described by the author. Ob- 
servations were then made of the turbulence activity by noting the direction of 


8 December, 1959, by Ronald T. Mc Laughlin, Jr. : 
_ 1 Research Engr., Iowa Inst. of Hydr. Research, and Asst. Prof., Dept. of Mechanics 


ind Hydr.; State Univ. of Iowa, lowa City, Iowa. : 
2 Research Assoc., Iowa Inst. of Hydr, Research, State Univ. of Iowa, Iowa City, Iowa. 
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motion of the fluid droplets. Three minutes after pouring, it was noted that 
about 30% to 40% of the fluid droplets were moving upward against gravity. 
After 7 min, less than 5% of the particles were moving upward, but the effect 
was still noticeable. Thus, there seems to be reasonable evidence that some 
turbulence effect would exist for about 10 min after pouring the sample. It 
would follow therefore, that some doubt might be raised as to the applicability 
of the author’s Eq. 17 


in computing the local mean settling velocity, because a diffusion term should 
be included in the equation. The appropriate equation would be, 


ag 
4 UN Gye aie 
09, 6 9) (te 5s) _, i. peer eee (2) 


ot oz az 


Integration of this equation with respect to z will yield, 
D 
a ao _= ao 
i Be OSes -Parge ee em ee a ie (3) 


instead of Eq. 18, as given by the author. A numerical solution of the integral 
may be obtained by determining the slope, for a particular depth, of a plot of 
the area under the concentration curve against time, but the numerical solution 
will not discriminate between the author’s Eq. 18 and Eq. 3. Thus, if the nu- 
merical result is divided by the concentration, an apparent fall velocity—here 
called Wa—will be determined which will approach the author’s w only in the 
absence of turbulence or as the turbulence decays sufficiently to cause the 
diffusion term to be insignificant. Assurance that the diffusion term of Eq. 3 
is important initially, may be obtained by considering that in order to obtain 
an initially uniform dispersion of sediment, the following equation must apply 


ag 
a(w 9) a(n ) ? 
Se ee Neve (4) 


wherein by definition the two terms must be equal. As the turbulence decays, 
the diffusion term will become less important and will cause Wa to rise to a 
peak and then decrease as a function of time in much the same manner as re- 
ported by the author (Fig. 8) for flocculation. It follows therefore, that the ef- 
fect of decaying turbulence on computations of the local mean settling velocity 
cannot be readily distinguished from the effect caused by flocculation withou! 
some other independent check. It is also interesting to note that the peak of the 
curve shown by the author (Fig. 8)occurs about 10 min to 12 min after the in- 
itial pouring which is the same order of magnitude as the time necessary fo 
the turbulence effect to become negligible, as estimated by the writers. 
Aside from the explicit effect of turbulence, an implicit effect may also oc: 
cur; that is, the settling of individual grams will be hindered—at a decreasing 
rate after pouring—as a resultof turbulence. Decaying, hindered settling woulc 
also tend to cause a peak to occur on a fall velocity - time graph. In order tc 
visualize these various effects, a schematic illustration, showing plots of fal 
velocity versus time for a particular sampling depth is given in Fig. A. 
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Each of the reported turbulence effects—decaying turbulence or decaying 
lindrance of settling—would tend to cause the lines of constant Z/t to slope 
cowards the Z axis of a concentration - depth graph (Fig. 7) and therefore be- 
nave similar to a suspension in which flocculation occurs. Because of this and 
he previously stated arguments, there is considerable doubt in the minds of 
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FIG. A.—EFFECT OF TURBULENCE, HINDERED SETTLING, AND FLOCCULATION 
ON COMPUTATIONS OF THE LOCAL MEAN SETTLING VELOCITY. 


he writers as to whether the author’s experiment isolates the effect of floccu- 
ation. Presumably, at great depths, turbulence would decay before an appre- 
iable change in the concentration gradient would occur, thereby reducing the 
mportance of including a diffusion term in the continuity equation. Perhaps 
his would enable the separation of flocculation from the other effects, but in 
der to establish the proper depths for sampling, further experimentation 
vould be required. 


_E. J. HALL.3—This paper is animportant addition to the literature on sedi- 
nentation and serves to narrow the gap between the theoretical approach to the 
rocess and practical design. 

“It is common to all the theoretical approaches from the time of Hazen that 
asic assumptions on the hydraulic flow in a tank and the settling character - 
stics of the suspension have to be made which are necessarily simplified to 
he point that practical application of the resultant theory is of doubtful value. 

The settling characteristics of the suspension, in particular, is a factor 
yhich has not received the attention whichit warrants and the common assump- 
ions made that the particles have a uniform settling rate, or that they settle 
rith uniform acceleration, is far removed from the truth, particularly where 
ewage is concerned. The graphs in this paper clearly illustrate this point. 


if Chf, Planning Engr., City Engr.’s Dept., Johannesburg, S. Africa. 
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The multiple depth pipette analysis suggested, deals adequately with this 
factor, but its application is likely to be difficult or at the least laborious in 
practice. 

In the search for a simple practical method of determining the settling 
characteristics of a sewage, the writer has made a close study of the process 
as it takes place in the laboratory and a brief report of the developments may 
be of interest. Observations show that even after a sewage sample has been 
standing for an hour or two, flocs can be seen settling at such high rates that 
they settle from near the surface to the bottom in a few minutes. 

It is presumed from these observations that settling of sewage proceeds as 
follows: 


1. Coarse material settles out with a few minutes quiescence only. 

2. The fine material remains in suspension slowly forming flocs. 

3. Once these flocs reach a critical size they commence to settle, and their 
removal in a shallow vessel is comparatively rapid. 


On this basis, the selection of an experimental method of determining the 
settling characteristics of a sewage, can be made on the assumptions that: 


a. The removal of the coarse material in sewage is easily achieved and 
may be ignored in tank design (except in as far as sludge storage and removal 
is concerned). 

b. The fine material remains in suspensionfor an appreciable period while 
flocculation proceeds. 

c. Once flocs of a critical size are produced, their removal, compared to 
the times required for their formation, is instantaneous. 


Based on these assumptions, therefore, samples taken from the surface of 
quiescent sewage at intervals as settlement proceeds will reflect the settling 
characteristics of that sewage. In the application of the method to practice, it 
may be added that as far as the instantaneous removal of flocs, once formed, 
is concerned, the assumption can be justified if the overflow rate of the tank 
under consideration is less than the settling velocity of the floc. 

The apparatus used by the writer to determine the settling characteristics 
of a sewage, consists of a drum 11 in. in diameter and 17 in. deep to which a 
small floating draw off, of 1/4 in. bore tube suspended on two corks, is fitted. 
A batch of sewage is added and samples removed from the surface at inter- 
vals. The samples are analysed either for suspended solids or preferably for 
P.V. (4 hr O.A.) as this is of more significance in works design. The experi- 
mentis repeated many times and the results averaged and represented graphi- 
cally as quality of effluent (or surface liquor) against time. The method ot 
analysis is well adapted to demonstrate the value of mechanical flocculation, 
as the batch after adding to a test vessel may be agitated in the vessel for < 
given time, at given paddle speeds, after which the paddle mechanism is re- 
moved, the floating draw off replaced and samples taken in the usual way. The 
use of two identical vessels simultaneously enables a direct comparison to be 
made between pure quiescence and flocculation followed by quiescence. Fig. 1 
shows a typical settling characteristics analysis curve of a domestic sewage 
and the characteristics of the same sewage when agitated by an oscillating 
paddle for 30 min at a maximum paddle velocity of 4.5 ips. 

It would appear, on first sight, easy to determine the efficiency of a sedi- 
mentation tank by relating its effluent quality to time of quiescence (or effec: 
tive detention period) by reference to the settling characteristics curve. Vari. 
ations in settling characteristics and sewage strength during the day are such 


TY 6 DISCUSSION 113 


lowever, that a reliable measure can only be obtained by averaging large num- 
vers of readings. A further factor is that flocculation does not proceed at a 
ixed rate but is greatly speeded up by slow mechanical agitation. Consequent- 
y, it must occur in varying degrees in all sedimentation tanks. For a graph 
¢ which direct reference could be made, it would, therefore, be necessary to 
sither: 


1. Determine the extent and duration of the agitation occurring in the tank 
nd duplicate this in the vessel used for analysis or 
2. Determine the optimum conditions for flocculation in the vessel and thus 
lbtain a curve which could be assumed torepresent the optimum and ideal con- 
ition. 
390 
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FIG, 1.—SETTLING CHARACTERISTIC ANALYSIS CURVE 


‘It is obvious that even with the simplified approach outlined previously, the 
rediction of tank performance from basic data bristles with difficulties in 
yhich unknown factors abound. It is to be hoped that further research such as 
hat described in this paper will be carriedout to bring more light to this fasci- 
ating problem. 
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TIDAL CHARACTERISTICS FROM HARMONIC CONSTANTS? 


Discussion by Remig A. Papp 


REMIG A. PAPP,! M. ASCE.—The paper fills a gap in our information about 
ides and the author should be congratulated for it. 

I would like to add supplementary comment to the section: Meteorological 
ffects. In addition to the winds, the barometric pressure can also influence 
1e height of the tide. After high barometric pressure the tide will be lower 
man computed and predicted; after a low pressure period the tide will be high- 
r. The difference could be quite appreciable and may amount to 6 in. to 8 in. 

This writer has participated in the port construction of the Port of Dun- 
uerke, France, in 1931-32. For the construction of the Jetee Est, 40 rein- 
xrced concrete caissons, precast in the inside of the harbor, had to be towed 
ut to the end of the jetty in construction. Due to the draft of the caissons this 
ywing could be made only twice a month during highest tide and every major 
eviation from the predicted tide level could have been critical. 

The field office has had a registering barometer and the differences between 
1e predicted andactual tide elevations were plotted under the barometric pres- 
ure curve. By plotting the plus deviation downward and the minus deviation 
pward, the two curves became similar. The tide curve, however, had had a 
ronounced lag behind the barometric curve which might amount to 10, 12, or 
ore hours. A difficulty is that the differences or deviations in the tide could 
e observed at the high and low points only since the other points are generally 
ot computed, having no importance. 

The notes of the observations were lost during the war but they were not 
yng enough to draw too many conclusions from them. It would be, however, 
uite easy to establish this type of observation in a few harbors with a high 
de range and to evaluate the influence of the barometric pressure on tide ele- 
ation. It is conceivable that after collecting enough data, a relation between 
1e angle of slope of the barometric pressure curve and the deviation in tide 
levation can be established for each local condition. 

In narrow estuaries the wind effect and the influence of the barometric pres- 
ure may have the same sign and may have to be superimposed or they may 
ave opposite signs. In such cases the direction and force of the wind should 
lso be taken into consideration. 


a December, 1959, by Bernard D. Zetler. 
dl Assoc., Hazen and Sawyer Engrs., New York. © 
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EARLY HISTORY OF HYDROMETRY IN THE UNITED STATES 
Discussions by William G. Hoyt and E. Shaw Cole 


WILLIAM G. HOLT, ! M. ASCE. —Mr. Kolupaila uses the beginning of World 
Var I (1914) as the end of the pioneer stage in the field of water measurement 
nown as stream gaging. In order to complete his listing of pioneers in this 
ield it seems appropriate to include other engineers who added to the know- 
sdge of the field in those early days. The writer refers to the following who 
ere district engineers of the United States Geological Survey, Dept. of the In- 
srior (USGS), in 1918 and who were not otherwise mentioned in the paper: 

C. E. Ellsworth, in Texas; A. L. Baldwin, in Idaho; F. F. Hemshaw, in Ore- 
on;O. W. Hartwell, in New England and New Jersey; A. H. Horton, in the Ohio 
‘iver Basin and who initiated the first nationwide system of collecting esti- 
1ates of developed and potential water power; William Lamb, in Montana; H. 
). McGlashan, in California; G. L. Parker, in Washington state and later chief 
ydraulic engineer (1939 to 1946); C. G. Paulsen, in the south Atlantic states 
nd later chief hydraulic engineer (1946 to 1957); A. B. Purton in the Great 
asin; R. C. Rice, in Kansas; and G. C. Stevens, in the Middle Atlantic region. 

Mr. Kolupaila suggests that those who measure water might well be called 
hydrometrists.” Perhaps the art of stream gaging has progressed tothe point 
here such a designation would be appropriate. The engineers referred to in 
1e paper or listed herein were far more than hydrometrists. Had they been 
nly what the definition implies, the Water Resources Division of the USGS 
ould not be the authoritative agency on national water resources it is today. 

Prior to World War I, the actual time spent measuring streams took but a 
mallpart of each 24 hr, a fact towhich the writer can personally attestas one 
f the pioneers listed by Mr. Kolupaila. When not traveling by train, livery, or 
n horseback, stream gagers were carpenters building structures; laborers 
hopping ice; public relation experts keeping farmers, ranchers, and others 
appy and reading gages twice a day for a maximum of $5 per month; comput- 
ig gage heights, preparing rating curves; making monthly reports to Washing- 
yn; and above all, exercising expert salesmanship by convincing state and city 
fficials, heads of irrigation districts and hydropower companies, and any 
thers who would listen, of the need for streamflow records. Without the fi- 
incial support thus obtained and the enthusiasm engendered, the work would 
xt have been maintained and certainly would not have continued at an ever- 
spanding rate to the present time. With but small staffs and large regions to 
ver, those in charge of the few district offices managed to maintain some 
000 gaging stations at a total cost of approximately $300,000, of which near- 
one-half was provided by state, public, and private agencies. 

In some quarters these early engineers have been criticized for not taking 
me to make a myriad of other observations deemed today as essential for 


2 January, 1960, by Steponas Kolupaila. a 
: Cons., Senate Select Committee on Natl. Water Resources, Washington, D. C. 
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proper water management and also for not thinking more profoundly about the 
hydrology of the river basins and the hydraulics and morphology of the streams. 
In general, it can be said that they were largely stream gagers, pure and sim- 
ple, but the records they recollected by sticking to their respective assignments 
were the beginning of a storehouse of streamflow information that now covers 
some. 200,000 station-yr of records, with annual expenditure by the USGS for 
all phases of water information now running about $25,000,000, of which co- 
operating bodies—local, state, and federal—still contribute about one-half. 

Unfortunately, at the beginning of the century there was no immediate need 
(and the future need was not adequately recognized) for information about ground 
water, the quality of water, the movement of sediment and other hydrologic data 
now. considered essential in the management of water resources. Even as ex- 
tensive as all the information on water resources is today, it is still far from 
sufficient to meet present-day needs if information now being collected by the 
Senate Select Committee on National Water Resources, of which Senator Robert 
S. Kerr of Oklahoma is chairman, is at all indicative. 

The writer is having an opportunity to analyze the extensive reports being 
received from the governors of the 50 states in response to a request from 
Senator Kerr for information as to present water problems and those antici- 
pated between now and 1980, together with recommendations as to what shoule 
be done toward their solution. Also, being reviewed is the testimony of some 
750 witnesses at 23 public hearings held throughout the country and numerous 
reports from federal agencies, consultants, and others. 

Although the writer is somewhat familiar with the availability of hydrologic 
data and deficiencies therein,2 he has been truly amazed with the apparent 
need disclosed in these sources for more hydrologic and meteorologic infor- 
mationof all kinds if solutions to the water problems of today and those antici- 
pated in 1980 are to be met. The extent to which this information is collected. 
analyzed, interpreted, and disseminated will have a profound effect on the ra- 
pidity and efficiency with which present and anticipated water problems will be 
solved. 


E. SHAW COLE,? F. ASCE. —There is some confusionin the account of thé 
pitot tube as used to measure the flow in pipes. A number of published paper: 
on the subject have sought to establish that there are two types, the simple ani 
the combined. The simple pitot tube is essentially a bent tube with the nozzl 
or orifice facing the flow. The total pressure against the orifice is a combi: 
nation of velocity headplus static headand to measure velocity in a pipe unde! 
pressure the static pressure must be measured by a wall piezometer. This i 
the type of tube used by Mills and Freeman. 

The combined type of pitot tube included a second tube that measured th 
static pressure. This was generally a long tube that was suited to open chan 
nel flow measurements. The Pitometer was a special type of combined tub 
that used two identical orifices 180° apart, one facing upstream and the othe’ 
downstream. This was developed and patented by Edward S. Cole and giventh 
coined name that included Henri Pitot’s name. Its special features were it 
compact size that permitted it to be installed easily in a pipe under pressure 
thus making it portable, its accuracy and its reversibility which provided - 
simple means for checking zero. 

2 «Water Facts For the Nation’s Future,” by Langbein and Hoyt, Ronald Press Co 
New York, 1959, | 

3 Pres., The Pitometer Assocs., New York, N. Y. 
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Two papers by power Ss. Cole, 4,5 one by C. W. Hubbard, § and the Power 
“est Code for 1949,‘ are important references for any one interested in the 
evelopment or use of the pitometer. 


7 
4 “Pitot Tube Practice,” by Edward S. Cole, Transactions, ASME, Vol. 57, 1935, p. 


31 
5 «pitot Tubes in Larger Pipes,” by Edward S, Cole and E. Shaw Cole, Transactions, 


SME, Vol. 61, 1939, p. 465. 
ba Investigation of Errors of Pitot Tubes,” by C. W. Hubbard, Transactions, ASME, 


st, 1939. 
“Power Test Code for Hydraulic Prime Movers,” 1949, pp. 135-165. _ 
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Journal of the Hydraulics Division 


Proceedings of the American Society of Civil Engineers 


february, 1960 
p. 41. The letter a at the end of Eq. 3 should be made a superscript. 


p. 41. In thetabulation near the bottom of the page the less than sign should 
e moved closer to the 1/2 in the first line. 


p. 47. In the equation following Eq. 11 the term 21 should be changed to 2}. 


p. 51. In line 11 the words “and the stage of the floodplain” should be 
hanged to “and the state of the flood plain.” 
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. PROCEEDINGS PAPERS 


> vhe technical papers published in the past year are identified by number below. Technical-division 
rshipis indicated by an abbreviation at the end of each Paper Number, the symbols referring to: Air 
Teasport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
es (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (SA), Soil 
shanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and Harbors 
, divisions. Papers sponsored by the Department of Conditions of Practice are identified by the sym- 
{PP). For titles and order coupons, refer to the appropriate issue of “Civil Engineering.” Beginning 
Volume 82 (January 1956) papers were published in Journals of the various Technical Divisions. To 
aeate papers in the Journals, the symbols after the paper number are followed by a numeral designating 
We issue of a particular Journal in which the paper appeared. For example, Paper 2270 is identified as 
279(ST9) which indicates that the paper is contained in the ninth issue of theJournal of the Structural 
WWision during 1959. 
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UNE: 2048(CP1), 2049(CP1), 2050(CP1}, 2051(CP1), 2052(CP1), 2053(CP1), 2054(CP1), 2055(CP1), 2056 
“HY6), 2057(HY6), 2058(HY6), 2059(1R2), 2060(IR2), 2061(PO3), 2062(SM3), 2063(SM3), 2064(SM3), 2065 
{ST6), 2066(WW2), 2067(WWa2), 2068(WWw2), 2069(Ww2), 2070(WW2), 2071(WW2), 2072(CP1)°, 2073(IR2)°, 
2074(P03)°, 2075(ST6)°, 2076(HY¥6)¢, 2077(SM3)°, 2078(ww2)°. 


ULY: 2079(HY7), 2080(HY7), 2081(HY7), 2082(HY7), 2083(HY7), 2084(HY7), 2085(HY7), 2086(SA4), 2087 
“(SA4), 2088(SA4), 2089(SA4), 2090(SA4), 2091(EM3), 2092(EM3), 2093(EM3), 2094(EM3), 2095(EM3), 2096 
‘(EM3), 2097(HY7)°, 2098(SA4)°, 2099(EM3)°, 2100(AT3), 2101(AT3), 2102{AT3), 2103(AT3),2104(AT3), 
2105(AT3), 2106(AT3), 2107(AT3), 2108(AT3), 2109(AT3), 2110(AT3), 2111(AT3), 2112(AT3),2113{AT3), 
“2114(AT3), 2115(AT3), 2116(AT3), 2117(AT3), 2118(AT3), 2119(AT3), 2120(AT3), 2121(AT3),2122(AT3), 
2123(ATS), 2124(AT3), 2125(AT3). 


\UGUST: 2126(HY8), 2127(HY8), 2128(HY8), 2129(HY8), 2130(PO4), 2131(PO4), 2132(PO4), 2133(PO4),2134 
4SM4), 2135(SM4), 2136(SM4), 2137(SM4), 2138(HY8)°, 2139(PO4)° 2140(sMa)°. 


EPTEMBER: 2141(CO2),2142(CO2),2143(CO2), 2144(HW3), 2145(HW3), 2146(HW3), 2147(HY9), 2148(HY9), 
2149(HY9), 2150(HY9), 2151(IR3), 2152(ST7)°, 2153(IR3), 2154(1R3), 2155(IR3), 2156(IR3), 2157(1R3), 2158 
(RS), 2159(IR3), 2160(IR3), 2161(SA5), 2162(SA5), 2163(STT), 2164(ST7), 2165(SU1), 2166(SU1), 2167(WWwa), 
“218e(ww3), 2169(WW3), 2170(WW3), 2171(Ww3), 2172(WW23), 2173(WW3), 2174(WW3), 2175(WW3), 2176 
“(wwa) 2171(WW3), 2178(CO2)°, 2179(1K3)°, 2180(HW3)®, 21¥1(SAS)°, 2182(HYS)°, 2183(SU1)®, 2184 
(ww3)°, 2185(PP2)°, 2186(ST7)°, 2187(PP2), 2188(PP2), 


CTOBER: 2189(AT4), 2190(AT4), 2191(AT4), 2192(AT4), 2193(AT4), 2194(@MA), 2195(EM4), 2196(EM4), 
97(EM4), 2198(EM4), 2199(EM4), 2200(HY10), 2201(HY10), 2202(HY1G), 2203(PL3), 2204(PL3), 2205 
PL3), 2266(PO5), 2207(P95), 2208(PO5), 2209(PO5), 2210(SM5), 2211(SM5), 2212(SM5), 2213(SM5), 2214 
), 2215(SM5), 2216(SM5), 2217(SM5), 2218(ST8), 2219(ST8), 2220(BM4), 2221(ST8), 2222(ST8), 2223 
, 2224(HY10), 2225(HY10), 2226(PO5), 2227(PO5), 2228(PC5), 2229(ST8), 2230(EM4), 2231(EM4), 
1232(A.T4)©, 2238(PL3)°, 2234( EM4)©, 2235(HY10)°, 2236(SM5)°, 2237(ST8)°, 2238(PO5)¢, 2239(ST8), 2240 

(PL3). 
OVEMBER: 2241(HY11), 2242(HY11), 2243(HY11), 2244(HY11), 2245(HY11), 2246(SA6), 2247(SA6), 2248 
SAG), 2249(SA6), 2250(SA6), 2251(SA6), 2252(SA6), 2253(SA6), 2254(SA6), 2255(SA6), 2256(ST9), 2257(ST9), 
1258(ST9), 2259(ST9), 2260(HY11), 2261(ST9)¢, 2262{ST9), 2263(HY11), 2264(ST9), 2265(HY11), 2266(SA6), 

(SA6), 2268(SA6), 2269(HY11), 22'70(ST9). 

ECEMBER: 2271(HY12)°, 2272(CP2), 2273(HW4), 2274(HW4), 2275(HW4), 2276(HW4), 2277(HW4), 2278 
fHW4), 2279(HW4), 2280(HW4), 2281(IR4), 2282(IR4), 2283(IR4), 2284(I1R4), 2285(POS), 2286(PO6), 2287 
06), 2288(P08), 2289(P06), 2290(PO6}, 2291(PO6), 2292(SM6), 2293(SM6), 2294(SM6), 2295(SM6), 2296 
M6), 2297(WW4), 2298(WW4), 2299(WW4), 2300(WW4), 2301(WW4), 2302(WW4), 2303(WW4), 2304(HW4), 
0), 2306(CP2), 2307(CP2), 2308(ST10), 2309(CP2), 2310(HY12), 2311(HY12), 2312(PO6), 2319(PO6), 
1A4(ST10), 2315(HY12), 2316(HY12), 2317(HY12), 2318(WW4), 2319(SM6), 2320(SM6), 2321(ST10), 2322 
O), 2323(HW4)°, 2324(CP2)¢, 2325(SMG)C, 2326(WW4)c, 2327(IR4)C, 2328(PO6)°, 2329(ST10)°, 2330 

vf 


VOLUME 86 (1960) =< << 
JARY: 2331(EM1), 2332(EM1), 2333(EM1), 2334(EM1), 2335(HY1), 2336(HY1), 2337(EM1), 2338(EMI 

39(HY1), Zantee ), 2341(SA1), 2342(2M1), 2343(SA1), 2344(ST1), 2345(ST1), 2346(ST1), cere 
B(EMIJC, 2349(HY 1°, 2350(ST1), 2351(ST1), 2852(SA1)¢, 2353(ST1)°, 2354(ST1). 

2356(CO1), 2357(CO1). 2358(CO1), 2359(CO1), 2360(CO1), 2361(PO1), 2362(HY2), 5 
, 2365(SU4). asia ser aU Boon eto pats AS 

(PC SM1), 2374(HY2), 2375(PO1), 2376(HY¥2), 2377(CO1)°, 379(SU1), 2380(SU1), 

( eed esealeUE y eeuare), 2385(SU1)©, 2386(SU1). 2387(SU1), 2388(SU1), 2369(SM1), 


JBRUARY: 2355(CO!1) 


| Sees Reece 
A 2(H : , ecient 2416(HW1), 2417(HW3), 2418(HWI)C, 2419WwI)e, 
1) ), 2424(SA2), 2425(8A2)°, 2426(HY3)¢, cypher e > or 
2432(P02), 4 , Pm 
2441(ST4), 2442(SM2), ry il pe 2 
A50(ST4), 2451(HY¥4), 2452(HY4), 
ey 2s08CP. 2466(ATI), 2467(AT1), 2468(SA3), 


| 2475(ST5), 2476(8A3), 2471(ST9), 
‘24B4(HYS), 24B5(HYS), 2486(AT1)C, 


